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INTRODUCTION 

This final report covers the period April 1, 1977 to December 31, 
1977 on Contract NAS9-14899-3S, 

Shown in Fig. 1 is a general flow diagram of work being done on this 
project. The work can be divided into two parts — an evapo transpiration 
model and a growth model. The end product of each is final grain yield 
of winter wheat,. The inputs into both models are identical — solar 
radiation, temperature, precipitation, axd leaf area index (LAI) . For the 
work considered here, the meteorological inputs are obtained from the 
National Weather Service but all have the potential for being estimated 
from spacecraft. Leaf area index (ratio of green leaf area to soil area) 
is obtained from ground and remote sensing estimates. Landsat computer 
compatible tapes (CCT) are used to estimate LAI. 

The first chapter of this report discusses the estimates of LAI from 
Landsat data. The yield estimates from evapo transpiration and growth 
models are discussed in chapter two and three, respectively. 
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CHAPTER ONE 

ESTIMATING. LEAF AREA, INDEX OF WHEAT 
USING LANDSAT DATA 
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ESTIMATING LEAF AREA INDEX OF WHEAT USING LANDSAT DATA 

Leaf area index (LAI) , the ratio of the plan area of leaves to a 
given area of land, is a very useful measurement. It is used in estimating 
evapo transpiration and yield (Kanemasu et al., 1977), and is an indication 
of plant vigor, disease, soil salinity and moisture stress (Knipling, 1967). 
However, it is an expensive measurement, requiring large inputs of labor 
and a large number of samples when done properly. Sampling errors can be 
quite large especially when sampling from nonuniform crop stands. 

With the launching of Landsat I in 1972 and Landsat II in 1975, a 
new' tool is available to the' agricultural community. Landsat I and. II 
have the capability to examine the. same area of land at the same local 
time every 18 days. The data relayed back to earth from Landsat consists 
of relative spectral intensity values from four wavebands from .4 pm to 
1.1pm (Table 1). The spatial resolution of Landsat is about .5 ha. Since 
most agricultural crop fields in Kansas are considerably larger than the 
spatial resolution, the spatial resolution is not a major problem; however, 
the temporal resolution of Landsat (18 days) is a concern in many agricultural 
studies. Because the MULTISPECTRAL SCANNER (MSS)' data represent average 
values for the entire field, sampling errors are reduced. The purpose of 
this paper is to examine Landsat-vegetation indices that are in the 
literature and to correlate them to ground observations of LAI. Then, we 
will simulate the LAI of wheat fields employing vegetation indices. 

Background . The atmosphere contains water, dust and particulate 
matter that vary in quantity from day to day and affect the quality and 
quantity of radiation detected by the MSS. A correction to the MSS measure- 
ments has been proposed for haze effects (MPAD, 1976) in which the smallest 
pixel value is subtracted from the other pixel values in an attempt to increase 
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Table. 1-. 

Spectral bands, for. the LAIJDSAT multi- 
spectral scanner (MSS) . 

Band 


Wavelength 

(cm) 

MSS 4 


0..5 to 0.6 

MSS' 5 


0.6 to 0.-7 

MSS- 6' 


O'. 7 to- 0 . 8 

MSS 7 


0 . 8 to 1.1 



contrast masked by the haze. Another type of correction for' haze and 
other materials in the atmosphere is to assinne all four bands are affected 
equally; therefore atmospheric effects are reduced by using ratios of the 
digital counts in each band rather than their absolute value. 

Variations in -the sun angle will also* affect MSS measurements. Some 
studies attempt to normalize the MSS values by dividing each band, by the 
sine of the solar elevation. Again this correction asstaaes that all MSS 
bands are affected similarly and band ratioing would correct sun angle 
effects. MPAD (1976) proposes that the sun angle has both additive and 
multiplicative affects for each band, the alpha values being the multipli- 
cative effects and the beta values the additive effects. A table- of alpha 
and beta values was developed for sun angles ranging' from 20 to 75 degrees 
for each individual band which allows the data to be transformed to a 
common sun angle. Sun angle affects the amount of radiation (lower amounts 
at lower sun angles) and the quality of radiation seen by the MSS. At lower 
sun angles, the radiation must penetrate more atmosphere to reach the 
surface; therefore, such effects as Raleigh scattering would be greater. 

Row orientation of crops in conjunction with differing sun angle causes 
variations in the amount of shading of soil and the lower parts of the crop 
(Fuch et al., 1972 and Suits, 1972); however, these effects are specific 
to individual fields and are difficult to take into consideration. 

In order to determine the vegetative properties such as LAI, the soil 
background must be taken into account. Wiegand and Gausman (1973) found 
that soils, regardless of color and moisture content, possess the same 
ratio between bands 5 and 7 and also between bands 5 and 6, Further, they 
show that on a diagram of band 5 versus band 7 (Fig, 1) , all foliage data lay 
on one side of the soil line while water, such as lakes and snow, lay on the 
other side of the line. In addition, the perpendicular distance (PVI) 




Fig-., 1. Diagram- of PVI (perpendicular vegeta-Cion index), 
showing 'the soil background line' and' the 
perpendicular line of a hypothetical point-. 
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between the soil line and the actual point under consideration is related 
to the foliage density. Therefore, soil effects could be effectively 
eliminated with the transformation, 

PVI = [((0.851*MSS5) + (0.355*MSS7))-MSS5)^ + 

((0.355*MSS5) + (0.148MSS7))-MSS7)^] [1] 
PVI6 = [(-0.498+(0.543*MSS5)+(0.498*MSS6)-MSS5)^ + 

(2.734+(0.498*MSS5)+(0.457*MSS6)-MSS6)^]'^ [2] 

where PVI and PVI6 are perpendicular vegetative indices; MSS 5, MSS 6, and 
MSS 7 are the digital counts in Landsat data. 

Kauth (1976) proposes another Landsat transformation in which he 
treats the four MSS bands as a vector. By transforming the vector with an 
experimentally obtained transformation matrix, the MSS data are converted 
into a four dimensional space — an axis in the greenness direction, soil 
brightness direction, yellowness direction and an axis in a direction he 
calls nonsuch (Fig. 2). The transformations can be expressed by a series 
of linear equations, 

GREENNESS (GVI) = 

-0.290 MSS4 - 0.562 MSS5 + 0.600 MSS6 + 0.491 MSS7 [3] 
YELLOWNESS (YVI) = 

-0.824 MSS4 + 0.533 MSS5 - 0.050 MSS6 + 0.185 MSS7 [4] 
SOIL BRIGHTNESS (SB) = 

0.433 MSS4 + 0.632 MSS5 + 0.586 MSS6 + 0.264 MSS7 [5] 
Kauth (1976) plots MSS data for corn fields on different soils in four- 
dimensional space and arrived with a figure he calls a "tasselled cap" due 
to its shape. This theory states that the greenness, soil brightness, and 
yellowness of a crop may be individually separated from MSS data. Thus, the 
greenness value (GVI) may be strongly related to foliage density alone- 
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Rouse et al. (1973) noted a correlation between (MSS7-MSS5) and green 
biomass. They developed an index called TVl (Transformed Vegetative Index), 
TVI = [((MSS7 -MSS5)/(MSS7 + MSS5)) + 0.5]'^ [6] 

and 

TVI6 = [((MSS6 - MSS5)/(MSS6 - MSS5)) + 0.5]’^ [7] 

In order to normalize (J1SS7-MSS5) they divided by (MSS5 + MSS7) , The 0.5 
was added to prevent the result from being negative and the square root 
was taken because of a better statistical correlation. 

Procedure . Since the launch of Landsat I in 1972, the Evapotranspiration 
Laboratory has obtained Landsat data and ground measurements of LAI on 
several winter wheat fields in Kansas (Kanemasu et al., 1974; Kanemasu et al., 
1977). Statistical correlations in MSS band ratios and LAI were performed 
(Kanemasu et al., 1977) and a regression equation was developed and given 
by 

LAI = 2.677-3.694(MSS 4/5) -2. 309 (MSS 4/6) + 5.751[MSS 4/(2x7)] + 0.43 

(MSS 5/6) - 2.692 [MSS 5/ (2x7)] + 3.071 [(MSS 4/5) - [MSS 4/ (2x7)] 
[MSS 4/5]. [8] 

In 1976, Kansas experienced, a very cool fall which resulted in reduced 
growth of winter wheat throughout the fall and winter and predictions of LAI 
from Landsat using [8] appeared to overestimate LAI. Presumably with a 
larger data set (115 observations) and the use of vegetation indices., a 
substantial improvement could be made and provide better results at low lAI 
where soil effects could predominate. 

Because of the large number of sources of variation involved (sun angle, 
atmospheric effects, LAI sampling, etc.) and the inability to quantify all 
of the possible sources, any practical model for estimating LAI must be 
empirical to some degree. Recognizing this, a linear model was developed 
using a statistical stepwise regression analysis.. The data used consisted 
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of 115 observations- of LAI and their respective Landsat MSS data. The 
various MSS band ratios and vegetation indices were individually correlated 
with LAI. The vegetation indices, showed about the same degree of correla- 
tion with LAI (Table 2), and MSS 5/6 was the most significant of the band 
ratios. The best-fit model was 

CLAI = -0.366 - 2.265 (MSS4/6) - 0. 431((MSS4/5) - (MSS4/7)) 

(MSS4/5)"^ + 1.7^5 (MSS4/5) h 0.057PVI [9] 

2 

where the R = .68.. It should be noted that MSS 4/5 appeared in [9] even 

though by itself it was poorly correlated to LAI (Table 2) . 

Shown iii Fig. 3 is the comparison of observed versus predicted LAI. 

At low LAI there was a large amount of variation presumably due to 

differences in soils, increased row effects, and other factors which 

predominate at lower densities. However, at LAI greater than about 0.5, a 

better correlation between LAI and. the model existed. Consequently, two 

models were developed — one for low LAI and one for high LAI. CLAI [9] 

was employed to classify which equation to use. The low LAI equation was 

developed by running the statistical program on all data with observed 

LAI less than 0.5. The high equation was obtained similarly with the high 

LAI data. The resulting equations can be expressed as: 

LAI = 1.093 - 1.138(MSS5/MSS6) - 0 .017 ( (MSS4/MSS5)-(MSS4/MSS7) ) (MSS4/MSS5) 

-0.016PVI6 CLAI £ 0.5 [10] 

and 

LAI = -5.33 + 0.036PVI + 6.54 TVI6 CLAI > 0.5 [11] 

Equation [9] incorrectly classified 7 data points as low, which should 

have been classified as high and did not misclassify any observed low LAIs as 

2 

high LAIs. Equation [10] performed poorly with an R of only 0.3 while [11] 

2 2 

had an R of 0.6. Overall the entire method had an R of .69 and is not 

2 

much better than [9] alone which had an R of .68. The justification for 
continuing to employ the 2-way classification system for predicting LAI is 
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Table 2 ... Correlation -of vegetation indices .and MSS' band, 
ratios with, leaf' area index (LAI) • 

2 

R- Values of MSS Transformations; 

2 2 
Vegetation Index R MSS Band,. Ratio R- 


PVI 

.55 

MSS 5/6 

-56 

PVI6 

.55 

MSS 5/7 

.48 

TVI 

.50 

MSS 4/7 

.27. 

TVI6 

.59- 

MSS 4/6- 

.20 

GVI 

.57 

MSS- 4/5 

.10 
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for the cases in which knowledge, of the ontogeny of the crop (crop calendar) 
allows one to determine the cases in which a field' is incorrectly classified. 
In these cases, the researcher can over-ride CLAI and the proper model can 
be employed. Although the procedure, is- empirical, the variability in both 
the satellite and the- ground measurements allows little chance for improve- 


ment . 
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CHAPTER TWO 

MODELLING DAILY DRY MATTER PRODUCTION AND 
YIELD OF WINTER WHEAT 



MODELLING DAILY DRY-MATTER PRODUCTION. AND YIELD OF WINTER WHEAT 


ABSTSACT 

Applxcability of many plant growth models are limited because of 
the input data requirements. Photosynthesis and respiration equations 
were developed from' meteorological data that could easily be obtained. 

The single crop parameter required was leaf area index. These equations 
were developed for winter wheat (Triticum aestivum , (L.)) from measurement 
of carbon dioxide exchange rate with field chambers and an infrared- 
gas-analysis system. Higher gross photosynthesis rates after jointing 
were attributed to sink enhancement of photosynthesis. Respiration was 
estimated as a photosynthesis-dependent growth component and a temperature- 
biomass dependent maintenance component. The equations predicted dry- 
matter production that agreed favorably with observed dry-matter accu- 
mulation. Head weight and yield equations were developed using predicted 
daily growth, LAI, and meteorological variables. Improvement of yield 
prediction was discussed. 

Key words ; 

Winter wheat, Gross photosynthesis. Respiration, Dormancy, Photo- 
synthetic efficiency. Sink capacity. Yield prediction. Grain to 
head weight ratio, CO^ exchange rate. 
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INTRODUCTION 

Plant growth and yield result from the interaction of environmental 
factors’, inherited traits of the plant, and current condition of the plant. 
Using models to study plant growth allows several factors to be varied at 
once. Models may be used to test theories of plant growth (Mllthorpe 
and Moorby, 1974), to estimate yield-, or to make crop production decisions. 
The simplest yield model involves summing or averaging several variables 
over the growing season and correlating those variables with yield by 
multiple regression. That approach requires a large number of crop cycles 
to develop a data base (Haun, 1974). More complex models may attempt to 
simulate individual plant processes on a daily basis (Milthorpe, Moorby 
and Morgan, 1974). For such models, equations to estimate daily photo- 
synthesis and respiration are essential. 

This study was designed to develop, for winter wheat ( Triticum 
aestivum (L.)), daily growth equations estimating photosynthesis and 
respiration based on carbon dioxide exchange rate (CER) measurements, 
and attempt to estimate yield from those equations for incorporation 
into an evapotranspiration-growth-yield model. Our objectives were: 

1, To investigate, whether or not photosynthesis and respiration 
change as functions of growth stage. 

2, To develop and 'test equations estimating photosynthesis, 
respiration, and dry-matter' accumulation as functions of leaf 
area index (lAl) and meteorological variables. 

3, To estimate yield from calculated CER values, and ontogenic 
and meteorological factors. 



METHODS' AND MATERIALS 


On 25 Sept. 1975, 1 ha plots o£ Triticum aestivum (L.) cvs. 

'Plainsman' and 'Centurk'- were planted, each at seeding rates of 34 kg ha 
and. 67 kg ha~^ at the Evapotranspiration Research Field 14 km southwest 
of Manhattan in Riley county, Kansas. On 1 Oct. 1974, 5 cultivars ('Centurk 
'Tamu', 'Prison', 'Sage', and 'Arthur') were planted at the site at 67 

I 

kg ha . On 22 Sept. 1976, eleven ^ ha plots were planted to varieties 
'Triumph', 'Eagle', and 'Plainsman'. Detailed measurements were made 
at the Evapotranspiration Research plots and less detailed measurements 
were made on commercial fields (in Riley, Ellsworth, and Finney counties) 
planted in late September to early October of 1974 and of 1975. 

From 25 Sept. 1975 to 24 Nov. 1975, and from 1 Mar. 1976 to 15 June 
1976 carbon dioxide exchange rate (CER) was measured with 4- open-chamber 
systems and two infrared-gas analyzers (liras II, Intertech Corp„). Every 
20 minutes each analyzer was switched between two chambers with an elec- 
tronically timed, solenoid valve system, previously described by Sij et 

2 

al. (1972). Each chamber enclosed an area of 1.2 m , including approx- 
imately 200 or 400 wheat plants, depending on the seeding rate; thus the 
CER values represented the effects of many plants. The chambers were 
moved to new locations at intervals (no longer than 14 days), depending 
on growth stage and weather. Soil was packed along the outside of each 

chamber to prevent large air leaks. The air circulation rate through the 

3 -1 

chambers, measured daily, was about 2.2 m min (air exchange about 1.2 
times per min). In an earlier study (Kanemasu and Hiebsch, 1975) dry 
matter estimated from chamber measurements agreed within about 10% of 
observed dry matter from jointing to heading. 

Based on measurements taken in the spring and fall of 1975, soil 



19 


respiration or carbon dioxide flux in mM CO^ m 12 hr from bare soil 
(without roots) was estimated as: 

Soil respiration - -16.4 - 2.61 (TMX -5- 2TMTN)/3 [1] 

during the night and 

Soil respiration = -34,34 - 1,48 (2 ti^B1X + PTMIN)/3 [2] 

during the day where TMX, THEN, and PTHEK are maximum and minimum daily 

temperatures and previous daily minimum temperature in C, These equations 

give values ranging between those found by Kanemasu et al. (1974-) during 

-2 -1 

summer at the Evapo transpiration Research field (-6.8 mM CO2 ra hr ) and by 

-2 -1 

Biscoe et al. (1975b) during, spring and summer in England (-2,5 mM 002® hr . ) . 

The number of chambers available precluded gathering replicated CER 
values'. The 20-minute readings were continuously recorded for each 
treatment for each day. Some data were discarded because of rapidly 
changing carbon dioxide concentration, probably caused by fluctuating 
atmospheric conditions (inversion with low wind). During May and June, 

1976, the solenoid valve from one chamber malfunctioned. 

To calculate leaf area index (lAI) and above ground dry matter 
accumulation we sampled plants weekly on the Riley county fields and 
approximately twice monthly on the Ellsworth and Finney county fields-- 
except during December, January, and February, when samples were taken 
less frequently on all fields. Leaf area was measured with an optical 
planiineter (Lambda Instruments). From heading to approximately hard- 
dough stage, 0,025 was added to lAI values to account for the photo- 
synthetic area of the heads (including awns). We measured yield by 
oven drying and weighing grain from three 6,1 m samples for each Riley 
county plot and one to three 3.05 m samples on all other fields except 
those in Finney county where U. S. Soil Conservation Service yield estimates 


Were used. 
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Solar, radiation (SR) was measured at Manhattan for the Riley county 
and Ellsworth county fields and at Dodge City for the Finney county fields. 
A thermograph was used to measure maximum and minimum temperatures at the 
Evapotranspiration Research field. At Garden City and Ellsworth, the 
National Weather Service temperature records were used. 

Light interception of the canopy was measured with one sensor (Lambda 
Instruments) pointing upward to measure incoming photosynthetic photon 
flux density (PPFD) , one sensor pointing down from' two meters height to 
measure reflected PPFD (RPPFD) , five sensors in parallel pointing upward 
at different positions underneath the canopy to measure transmitted PPFD 
(TPPFD) , and five sensors in parallel facing downward beneath the canopy 
to measure PPFD reflected from the soil (RSPPFD) . Intercepted PPFD (IPPFD-) 
was calculated as 

IPPFD = PPFD - RPPFD - TPPFD + RSPPFD ' [3’] 

-2 -1 

where PPFD is in micro Einsteins cm day . One Einstein is one mole 

of photons. For most of the growing season, RSPPFD was small. 

On fields where ontogeny data were not recorded (Ellsworth and Finney 

counties), stage of growth' (BMTS) was estimated using a biometeorological 

3/ 

time scale (Feyerherm and- Paulsen, 19-76)— (Table 1) where emergence has 
a value of' 1, jointing of 2, heading of 3, soft dough of 4, and maturity 
of 5. 

Because of the similarity of the pathways of photosynthesis and 
transpiration, the responses to water stress are similar. An evapotran- 
spifation model for sorghum and soybean (Kanemasu, Stone, and Powers, 

1976), adapted to winter wheat, was used to estimate the effect of water 
stress on photosynthesis. Gross photosynthesis (GP) was assumed to be 
reduced by a water stress factor (K^) : 

3/ • 

— Feyerherm, A. M. and G. M. Paulsen. 1976. A biometeorological time 

scale applied to winter wheat. Agron. Abstr. p. 9. 



Table 1. Equations used to calculate biometeorological time scale growth stages for winter wheat, 

with daylength in hours and maximum temperature (Tmax) and minimum temperature (piin) in F. 

VI = C(k,l) Daylength - c(lc,0) + C(k,2) Daylength - C(k,0)^ 

V2 = C(k,4) Tmax - C(k,3) + C(k,5) Tmax - C(k,3)^ 

V3 = C(k,6) Tmin - C(k,3) + C(k,7) Tmin - C(k,3)^ 

TDAY = VI (V2 + V3) 

BMTS = STDAY = SV1(V2 + V3) (BMTS < 2.00) 

BMTS remains less than 2.0 until TDAY > 0.02 10 times after BMi’S > 1.95 or until 2TDAY > 2.8. 

When BMTS == 2.0, then T^TDAY is set to 2.0. 


k 

Growth stage 

1 

2 

3 

C Matrix 
4 5 

6 

7 

8 

0 

Planting to emergence 

lo'^ 

-.1419 X 10"^^ 0 

44.37 

.07652 

-.001571 

.06875 

-.001597 

1 

Emergence to jointing 

8.413 

.05581 

0 

23.64 

-.06324 

.000905 

.006601 

-.771 X 10“'^ 

2 

Jointing to booting 

10.93 

.02613 

-.001701 

42.65 

.01047 

0 

.01396 

0 

3 

Booting to heading 

10.93 

.02613 

-.001701 

42.65 

.01047 

0 

.01396 

P 

4 

Heading to soft dough 

10.94 

.02021 

-.001192 

42.18 

. .01688 

0 

.02136 

0 

5 

Soft ripe to maturity 

24.38 

-.02165 

0 

37.67 

.003543 

0 

,01811 

0 



2'2 


K = 6/.350max 
s 

K = 1 
s 


(9<.358inax) 

C0>,359max) 


where 0 is, available soil water and 0 max is the available soil water at 
field capacity. 


RESULTS AND DISCUSSION' 


Model Development 


Daily dry matter production (DMP) can be estimated from the 24 hour 
CER or TGERo The CER measurements were corrected for soil respiration 
according to [1] and [2] and integrated to obtain the 24 hour CER estimate 
(TCER), TCER is then the algebraic sum of the daytime CER (DCER). and the 
negative nighttime CER, 

McCree (1974) estimated respiration for sorghum- and white clover with 
equations of the form 

Resp = a dCER + b dmA [6] 

where a is a growth coefficient and b is a temperature dependent main- 
tenance coefficient. We estimated the nighttime CER or dark respiration 
(NResp) as : 

NResp = (N^)(-,627 - 0,0148 DNCE - x DMA) [ 7 ] 

-2 -1 

where NResp is the nighttime flux of CO 2 (mM CO 2 ^ day ) ; N^ is the night 

length in hours'; DCER is gross photosynthesis minus the daytime' respiration 
-2 "1 

(mM CO 2 m day ) ; DMA is the- accumulated above- and below-ground dry 
2 

matter (gm/m ) until maximum LAI (MLAI) is attained; thereafter DM is 

multiplied by lAl/MLAl, Thus, when all leaves have senescenced, the 

maintenance component is zero. The "maintenance” coefficient (M ) is 

n 

temperature dependent and of the form 

M = 10 "^ X (529.3 + 1,857 x T + ,0529 x T ^) 
n n n 


[ 8 ] 
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where T is (T + 2T . )/3, T and T . are the 24 hour maximum and 
n max mrn max min 

minimum temperatures in degrees Kelvin. Equation [7] --evaluated by 

using chamber measurements of DCER, night respiration, obseirved dry matter, 

2 

and temperature data — resulted in an r = 0.7^. The constants in [7] and 

[8] were determined with a multiple regression program. 

Daytime respiration (DResp) can be expressed by [7] with daylength 

substituted for nightlength,. gross photosynthesis (GP) substituted for 

DCER and M, substituted for M . M, is given by [8 ] with T, = (2 t 

d n d d max 

+ FT . )/3 substituted for T ♦ PT . is the previous daily minimum temper- 
min n min 

ature in degrees Kelvin. 

Table 2 compares the maintenance and growth coefficients for wheat 
with those of McCree's over a range of temperatures. The wheat coeffi- 
cients were less sensitive to temperature than those of sorghum and white 
clover, probably because wheat is a cool season crop. 

Gross photosynthesis (GP) depends to a large extent on the number 
of photons or quanta of light intercepted by the canopy. Shown in Fig. 1 
is the relationship between intercepted PPFD (IPPFD) and leaf area index 
(LAI), IPPFD can be expressed as 

IPPFD = .5739PPFD LAI ^.6 [9] 

2 

with r = .87 and 

IPPFD = .95 PPFD ■ LAI>4o6 [10] 

-2 -1 

where IPPFD and PPFD are given in pE cm day . Estimating light inter- 
ception by the point quadrants method (L7ilson, 1967) while probably more 
accurate, would require measurement of leaf angle distribution or mean 
leaf angle, so was not used. 

When PPFD was not measured, PPFD was estimated from daily solar 
radiation (SR) by the relationship 


PPFD = 9.07 SR 


[ 11 ] 
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af 


Table. 2',. itei-ntenance -andl growth- coefficients compared, for sorghum' and • 

white, - clover. (=McCree, ..1974.)’ and. winter' wheat (from'.-0'°'' C to 30® C 
when daylengthr fs- 12.. hours). 


TEMPERATURE 

■ MAINTENANCE- COE-FF-IGI-ENTS 

- 

(c> 

Sorghum- 

White Clover. 

Wheat 

0 

' . 001 

.002' 

,0050. 

10; 

.003 

,008 

.0053 

20 

' .009 

.023 

, 005 6-- 

30- 

..018' 

,049 

,0059- 


GROWTH: COEEF-ICIENTS 
Sorghum- Wh-ita Clover: ’ Wheat 


0,f4 


0,14- 


0,18. 
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-1 

where SR is in ly day . Equation [11] was obtained from SR and PPFD 
measurements at the Evapo transpiration Research site. Therefore, from 
measurements of solar radiation and LAI, intercepted PPFD can be esti- 
mated. 

Gross photosynthesis (GP) is reduced by extremely hot or cold temper- 
ature. Based on our CER data photosynthesis seemed to be reduced by max- 
imum’ temperatures above 23C and by minimum temperatures below 5C. For hot 
days, photosynthesis is decreased as increases and air temperature 

increases above the critical temperature for an increased portion of the 
day. Also when the minimum' temperature is increased, for- a given maximum 
temperature, then the portion of the day with stress increases. A high 
temperature stress factor is calculated as : 


Tjj = [23 - 

1/2(T . + PT . )]/[T - 1/2(T . + PT . )] 

min min max min min' 

[12] 

when T > 23 C 
max, 

5C -< 1/2 (T . + PT . )< 23C 

min min — 



(T < 23C) 

max — 

[13a] 

Th = o 

[1/2'(T . + PT . ) > 23C] 

mxn mm 

[13b] 


Similarly low temperature stress increases with decreased maximum 
and minimum temperatures. 


T = (T - 5)/(T - PT , ) 

L max max mxn 




T = 0 

Li 


(PT . > 5°C) 

min ~ 

(T < 5°C) 
max — 


■OE SOOB 


q'UALX'5^ 


Temperature stress is the product of T and T : 

H L 


T = T T 
S H L 


[14a] 

[14b] 

[14c] 


[15] 


Potential GP (without water or temperature stress) was estimated 
for various growth stages by using the chamber and light-interception 
data. Following Biscoe et al. (1975a) monomolecular equations were used 
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to estimate the photosynthetic response to intercepted light. From 
emergence to jointing 

GPp = 1000 [1 - exp (-.000276 IPPFD) ] [16] 

2 

with r = .84, where GP^ is potential gross photosynthesis in mM CO^ 

“ 2 1 

m day . During cold days (T < 5 C) , GP =0 (Martin and Leonard , 

max p 

1967). From jointing to maturity and LAI ^1.5, 

GP^ = 2300 [1 - exp (-.00057 IPPFD)] 

2 

with r = .68. If LAI < 1.5, [16] is used. 

[17] were evaluated by regression analysis. 

is calculated as: 

GP = GP K T [18] 

p s s 

Table 3 illustrates the difference in photosynthesis, respiration, 

and daily dry matter production (DMP) (calculated by [20]) from the 

prejointing and postjointing photosynthesis equations. The difference 

in gross photosynthesis between equations [16] and [17] may be due to 

sink enhancement of photosynthesis (Evans and Dunstone, 1970) after 

jointing. GP was not allowed to be negative. 

P 

Estimated daily gross photosynthesis and respiration are summed 

in Table 4 by growth stage. Denmead (1976) estimated that a spring 

-2 

wheat crop had 31425 mM CO^ m daytime net photosynthesis and 12600 
-2 

mM CO^ m nighttime respiration from one week before heading to 26 

days after heading. For Plainsman, the values from heading to 35 days 

-2 

after heading were (Table 4) 9662 mM CO^ m daytime net photosynthesis 

-2 

and -2929 mM CO^ m nighttime respiration. Biscoe et al. (1975a) 

-2 

reported 23750 mM CO^ m daytime net photosynthesis and 5795 mM 
-2 

CO2 m nighttime respiration for spring barley for a similar period. 


[17] 

The constants in [16] and 
Actual gross photosynthesis 
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Table. 3, Model estimates of gross photosynthesis (GP), day and night 

respiration (Resp) , and daily dry. matter production (DMP) for 
a range of intercepted light values when. accumulated' dry 

2 * 
matter is 500 gm/m , ,Tm is 20 G, Tn. is 15 C and daylength is 12 hour: 




Prejointihg, 





II>AR 


n ■ 

Re sp- 

DMP — 

fi-E 

.-2 -1 
cm. day 

■ 

mM CO 2 

-2 , -1 
m day 

gm/m^ 


0 

1/ 

(000)-^^ 

0 

- 84 

- 2.4 


1814 

(200,) 

394 

-205 

5.8 


3628 

(400) 

633 

-282 

10. ,6 


4535 

(500) 

714 ‘ 

-308 

12.3 



Postjointing 





IPAR 


Gpl/ 

Resp-- 

DMP^"^ 

p-E' 

-2 --1 
cin day 


mM CO 2 

-2.-1 • 
m day 

gm/m^ 


0 

(000)-^ 

0 

- 84 

- 2.4 


1814 

(200) 

1482 

-557 

27.8 


3628 

(400) 

2009 

-111 

38.5 


5442 

(600). 

2197 

-788- 

42.4 


•^'^Intercepted light in ly day 

— G? = 1000'[1 - exp (-,.000276 IPAR) ] . 

•^'^GP = 2300;[1 - exp (-.00057 IPAR)]. 

4/ 

— Resp' =-NResp' + DResp. 

A/pMP = *03 (GP + DResp -t-NResp) = ,03(GP + Resp), 
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Table 4, Predicted gross photosynthesis (GP), daytime respiration 
(DResp), nighttime respiration (NResp) , and net photo- 

-2 -1 

synthesis (TCER-) in units of mM CO,,, m day and dry 

^ 2 

matter' accumulation (DM) in units of gm/m for Plainsman 
from emergence (1)- to jointing (2) , from jointing to 
heading (3), and from heading to soft dough (4), and number 
of days for each stage summarized. 



1-2 

GROWTH 

2-3 

STAGE 

3-4 

TOTAL 

Days 

189 

36 

36 

260 

GP 

22331 

57594 

12022 

91947 

DResp 

- 4444 

-13972 

- 2360 

-20776 

NResp 

- 7409 

- 8268 

- 2929 

-18606 

TCER 

10478 

35354 

6733 

52565 

DMA 

339 

1145 

.218 

1703 



For the season (emergence to soft dough), 42.8% of the gross photo- 
synthesis- was respired which is less than estimated by Biscoe et a-1. 
(1975b). and Connor (1975) of 48.6% and 55%, respectively. 

Photosynthetic efficiency values, in terms of total solar radiation 
(SR in ly/day), were, caiculated for Plainsman using the heat of combus- 
tion of plant material as 4Kcal/gm suggested by Lieth (1968.) as an 
appropriate average for herbaceous plants,, For Plainsman, values 
were 0.46%, 1.61%, and 0,79% respectively from emergence- to jointing, 
jointing to heading, and heading to soft ripe. From data of .Biscoe et al. 
(1975b) values for spring barley were 1.55% for the six weeks before 
anthesis and 1,00% from anthesis to soft ripe. All values were consid- 
erably below the estimate by Loomis and Williams (1963) of 5-6% for 
maximum possible efficiency. 

The 24-hour CER or TCER can be obtained from [7]-, [16], [17], 
and [18] which require measuring solar radiation, maximum temperature, 
minimum' temperature, and leaf area index, 

TCER = GP + DResp + NResp, [19] 

From TCER estimates, daily dry-matter production (BMP) with 8% ash 
can be calculated by converting mM CO 2 £0 gm of carbohydrate 

BMP = TCER X 30/1000 x 1.08 [20] 

2 

where BMP is in gm/m . By summing the BMP for each day, the cumulative 
value can be compared with total dry matter observed in the field. 

Tests 

To compare the values predicted by the model with measured above- 
ground dry matter, we made several assumptions about root dry-matter 
production. Root dry weights were taken in the 1974-75 season on 
Centurk (Kanemasu et al,, 1977). Rooting densities of llO, 180, and 



250 gm/m were obtained for dormancy, jointing and heading, respectively, 

Pasternak (1974), Cok and Wright (1975), and Connor (.1975) found rooting 

2 ‘ 

densities of 230, 220, and. 210 gm/m , respectively, at maturity. Those 
root values, added to the above-ground dry weights, were compared with 
those of the model for Riley county fields. Because seeding rates were 
lower in western than in eastern Kansas, we arbitrarily used 3/4 and 
1/2 of root weights for Ellsworth and Finney counties, respectively. 
Fortunately, an error in estimating root weights did not result in 
large errors' in total dry matter, especially late in the growing season. 

Shown in Fig, 2 are the model-predicted dry matter and observed dry 
matter for all fields in which data were not used in developing the model. 
In general, the values' lie along^ the 1:1 line. 

Shown in Figs, 3, 4, and 5 are the temporal trends in observed and 
predicted dry matter and leaf-area index. Shown in Fig, 3 are predicted 
and' observed dry matter for Evapotranspiration Research fields. Pre- 
dicted and' observed dry matter agreed closely until heading (1 May), 
when the equations began to underestimate dry -matter accumulation. The 
NCE measurements also decreased after heading as lAI decreased. That 
discrepancy, which did not occur on the commercial fields studied 
.(Figures 2, 4, 5), might have resulted from errors (1) in sampling dry 
matter; (2) in measuring NCE after jointing; (3) errors associated with 
photosynthesis by yellowing leaves (not included in LAI measurements); 
or (4) in estimating respiration after heading. After a severe- frost 
(3 Ifey), the chambers might have been moved to locations with greater 
head damage than the dry matter samples. 

No CER measurements were made during December, January, or February; 
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. 3, Seasonal trends in observed and predicted dry matter 
accumulation and LAI. for Evapotranspiration Research 

fields (1975-1976), 
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Fig. 4. Seasonal trends' in observed and predicted dry matter accumulation' 
and LAI for a Riley county field' (1974-1975) ,, 
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Fig; 5. Seasonal trends in observed and predicted dry matter 

accumulation and LAI for an Ellsworth county field (1975-1976), 




hence, we could not evaluate directly the assumption that photosynthesis 
ceases when the maximum temperature is less than 5 C, However, the 
predicted dry matter at the end of February was not greatly different 
from the observed dry matter (Figs. 2, 3, 4, 5). Repka and Kubova (1971) 
reported that roots and shoots of V7heat respond rapidly to temporary 
increases in temperature. 

We could not adequately test the assumption that GP decreases 
linearly with available soil moisture after a threshold value because 
water deficits did not occur until the ripening period « At that time, 

GP was also declining because of leaf senescence. 

Yield Process 

Carbohydrate enters the grain from carbon fixed in the leaves or 
head during the grain-filling period or from carbon fixed earlier and 
released for translocation to the head during senescense (Lupton, 1972). 
Although yield. is influenced by the' number of culms per area of field, 
determined during fall and spring tillering, yield formation actually 
begins with head formation at ear. initiation (jointing), when the 
potential number of seed-bearing tillers and ntimber of seeds per head 
are set. 

Five types of factors may be regarded as limiting yield, of which 
three (2, 3 and 4) are considered in the model, 

1) Sink capacity is set by tillering, panicle initiation, and 
flowering. Very high or low tiller number per area may be 
partly compensated by the percent of tillers forming heads and 
by the number of seeds per head. Stress during tillering, 
jointing, or flowering may reduce sink (Halse and Weir, 1974; 

2 

Gifford, 1977), -Sink capacity, then, may be defined as heads/m 
X seeds/head x potential weight/seed. 
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2) Sources contributing to grain filling are TCER after jointing 
and translocation of carbon, TCER is calculated as a function 
of intercepted light, temperature, growth stage, water stress, 

and temperature stress'. Carbon' released is estimated from decrease 
in lAl (reflecting senescense of leaves, sheaths, stems, and roots). 

3) Translocation rate of carbon from, leaves and stems into the 
grain may be limited by high temperature stress (Asana and Basu, 
1963). 

4) Duration of the grain-filling period is, estimated in the model 

from daylength, light intensity, and temperature by a biometeor- 

4/ 

ological timescale— (Table 1) but varietal differences are not 
considered. 

5) Damage to the crop during grain- filling due to biological or 
mechanical factors (pests, hail, flooding) is not considered, 

TThen sink capacity is not limiting, then rate of supply of the sources, 
rate of movement' into the grain, and duration of the grain-filling period 
determine yield. 

Yield Equations 

Head weight increases slowly from jointing until shortly after 
anthesis ; thereafter it increases rapidly in weight as the grain is 
filled. Summarized in Table 5 are changes in head weight and lAl during 
several periods from jointing to maturity, and TCER during each period, 
for 11 plots at the Evapotranspiration Laboratory Research site in 1977, 

,*47 

Feyerherm, A. M. and G. M. Paulsen. 1976. A biometerological time 
scale applied to winter wheat. Agron. Abstr. p. 9. 



Table 5. Summary of changes in head weight, total weight, , and LAI for the Riley county fields used to 
develop yield equations and coefficients derived for the yield equations. 



Stage 

2.0 

- 2.7-^ 


2.7 - 3 

.oi/ 


3.0 - 5.0^^ 



Variety 

Plot 


Head 

Growth 

Total 

Growth 

II / 
ALAI- 

Head 

Growth 

Total 

Growth 

:alai-'^ 

Head ^ . Total^ , 

Growth- Growth—' 

7/ 

DMA,- 

4 



(gm/m^) 

(gm/m^) 

(gm/m^) 

2 

(gm/m ) 

Triumph 

ET Lab 

Tl 

20 

136 

.19 

30 

18 

i .41 

292 

90 

387 

Triumph 

ET Lab 

T2 

20 

273 

.34 

20 

19 

: .84 

275 

104 

538 

Triumph 

ET Lab 

T3 

20 

183 

.33 

30 

19 

.67 

327 

84 

433 

Eagle 

ET Lab 

El 

20 

284 

.28 

30 

19 

.92 

361 

98 

650 

Eagle 

ET Lab 

E2 

20 

350 

1.14 

70 

96 

.54 

258 

86 

755 

Eagle 

ET Lab 

E3 

20 

257 

.35 

20 

19 

.70 

327 

93 

617 

Eagle 

ET Lab 

E4 

20 

173 

.15 

30 

14 

.16 

292 

70 

506 

Plainsman 

ET Lab 

P2 

30 

82 

.35 

30 

18 

.38 

275 

71 

385 

Plainsman 

ET Lab 

P3 

50 

83 

.42 

20 

18 

.39 

327 

69 

385 

Plainsman 

ET Lab 

P4 

40 

83 

.32 

30 

18 

.44 

292 

74 

390 

n.a. 

Puett 


20 

415 

.31 

50 

87 

1.36 

293 

136 

763 


2 : 


280 

2319 

3.54 

360 

345 

6.45 

3319 

975 

4980 

Coefficient 


.13 



.71 



1.00 


.41 


1 / 

— Jointing to booting stage. 

2 / 

— Booting to heading stage. 

3/ 

~ Heading stage to physiological maturity. 

4/ 

— AIAI values are limited to < 1,00 in any period for calculating head growth coefficients, 

5/ — 

~ Divided by 31 ^^ (average temperature stress from heading to physiological maturity) to correct for 


effects of temperature stress. 

6 / 2 

—Total growth is limited to < 300 gm/m for calculating coefficient. 

7/ 2 

— DMA at 4,0 is limited to < 500 gm/m for calculating coefficient. 


03 
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Assuming sink capacity is not limiting, head weight is calculated as: 


2 7 

" 30 

Head Weight = 1.08[ .13( S TCEfcj^) + 


3 3 

.71 ( A lAl X SLW + E TCER- 
2.7 2.7 


30 

1000 


)+ 


_ 5 5 

(^Ts^K A LAI X SLW + E TCER^^ + .ISDMA^)]- [21] 

2 ^30 

with r = .72, where 1,08 adds 8% ash; E TCE R^qqq is total CO'^ exchange 

a 

b 


30 


rate of the crop, from stage a to stage b, converted to carbohydrate; A LAI x 

a 

SLW is decrease in LAI from stage a to stage b times the specific leaf 
2 — 

weight (45 gm/m ; Puckridge and Donald, 1967); average temperature 

stress from stage 3 to stage 5; and DMA^ is total dry matter accumulation 
b 

up to stage 4, A LAI x SLW is used to approximate carbon released by 
a 

senescence of leaf sheaths, stems, and roots as well as leaves. To prevent 

overestimation of yield (possibly when sink is limiting) change in LAI was 

2 2 2 2 
limited to .5 m /m from stages 2.7 to 3.0 and to 1.0 m /m from stages 

30 2 

3.0 to 5.0, TCER-Jq^ was limited to 150 gm/m from stages 2,7 to 3.0 and to 

2 2 
300 gm/m from stages 3.0 to 5.0, and was limited to 500 gm/m . 

Temperature stress during grain-filling may affect yield by reducing 

sugar translocation from the stem into the head as well as by reducing 

the photosynthesis rate (Asana and Basu, 1963). Ruckenbauer (1975) 

suggested that translocation rate into the grain is dependent upon the 

capacity of grain to store carbohydrate-. So temperature may limit grain 

growth by reducing the starch synthesis rate in the grain. 

Grain and head weights for winter wheat and spring wheat from several 

sources are compared in Tig. 6 and in Table 6 with references. Most of 

the data from the literature lie close to the .74 grain/head ratio. Some 

of the points with lower ratios can definitely be attributed to late 
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Table- 6a. Comparison-, of grain- and head- weights- for s.e,v.eral -spring 
and -winter- wheat fields reported, in. the- Literature'. 


Reference 

Species 

Grain - 
Wt„ 

(gm/m ). 

Head 

^2-, 
(gm/m- ) 

■ Ratio 

Puckri'dge and- Donald', 196-7' 

Spring 

46 

65,2 

.71 



173 

236 - 

.73 



247 

342' 

.72 



234- 

2-99 

« 

00 



185 

249 

.76 • 

Spiertz',. 1974- 

Spring- 

80.1. 

1124 

.71 

Pearman-, Thomas,, and. Thorne, 1977 

Spring 

481 

6'50 

.74 



514 

685... 

.75 



444-. 

600. 

.74 

Paltridge- et'.al,. , 1-97-5, 

Winter . 

281.6- 

.360.9“ 

.78- 

Pal'tridge. et al. , 1.972^ 

Winter 

4o,r.o 

627,1 

.64- 
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Table 6b. 

Comparison of 
winter wheat : 

grain and head 
fields . 

weights 

for several 

Kansas 


County 

Field 

Variety 

Year 

Grain 

Wt 

Head 

Wt 

Ratio 

2 

(gtn/m ) 

Riley 

Lysiraeter Low 

Plainsman 

1-976 

195 

660 

.30 


ET Lab Low 

Centurk 

1976 

306 

700 

.44 


ET Lab Med- 

Centurk 

1976 

339 

680 

.50 


ET Lab High 

Centurk 

1976 

289 

680 

.43* 


Cedarburg 

n.a. 

1976 

183 

470 

.39 


Puett 

n.a. 

1976 

174 

480 

.36 


ET Lab Tl 

Triumph 

1977 

165 

210 

.78 


ET Lab T2 

Triumph 

1977 

159 

220 

.71 


ET Lab T3, 

Triumph 

1977 

183 

200 

.91 


ET Lab T4 

Triumph 

1977 

110 

190 

.57 


ET Lab El 

Eagle 

1977 

208 

260 

.80 


ET Lab E2 

Eagle 

1977 

183 

240 

,76 


ET Lab E3 

Eagle 

1977 

183 

2:10 

.88 


ET Lab E4 

Eagle 

1977 

153 

210 

.74 


ET Lab Pi 

Plainsman 

1977 

128 

200 

,65 


ET Lab P2 

Plainsman 

1977 

147 

220 

.67 


ET Lab P3 

Plainsman 

1977 

183 

250 

.73 


ET Lab P4 

Plainsman 

1977 

165 

220 

,76 


Cedarburg 

n.a. 

1977 

242 

400 

.61 


Puett 

n.al 

1977 

164 

250' 

.66 

Finney 

18 

n.a. 

1977 

206 

419 
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n.a. 

1977 

167 

225 

.74 


50E 

n.a. 

1977 

153 

440 

,35 


94 

n « 3- 0 

1977 

162 

290 

.56 


129 

n.a. 

1977 

185 

481 

.38 


186 

n.a. 

1977 

240 

360 

.67 
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frosts (Lysimeter, 1976), flooding during grain- fill (Plainsman 1,, 1977), 
and fertilizer burn (Triumph 4, 1977). 

Predicted yields, calculated as 

Yield = .74 X' Head Weight [22] 

2 

with r = .41 are plotted against observed yields in Fig. 7 and listed 
with observed yields and obseirved and predicted' head weights in Table 7, 

To improve the method of yield prediction outlined in this paper, 
we need information in several areas. Tillering occurs in response to 
light, temperature, population density, nutrients, and previous growth 
of the crop with different varieties responding differently to these 
conditions. Spikelet initiation and possibly seed setting are' affected 
by temperature stress (Halse and Weir, 1974), shading (Pepper and Prine, 
1972), and reduced leaf area (Davidson, 1965), at jointing and at anthesis. 
The greater variation in the grain/head weight ratio of wheat crops in 
Kansas than that of other wheat crops reported in the literature (Table 
6) may be due to high temperatures in Kansas in early summer. The 
different varieties planted in 1975, 1976, and 1977 in Riley county had 
jointing and heading dates different by as much as 6 days (Table A27 in 
Appendix A), differences not reflected in the- biometeorological timescale. 

CONCLUSIONS 

The maintenance coefficient of wheat presented here was less sensitive 
to temperature than were coefficients found for sorghum and white clover 
(McCree, 1974), possibly because wheat is a cool-season grass. Photosyn- 
thesis was found to increase sharply after jointing possibly because of 
increased sink capacity. On most wheat fields studied, predicted and ob- 
served dry matter agreed for most of the growing season. Assuming that 
photosynthesis ceased when the maximum temperature was less than 5°C gave 
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Table 7. Listing of observed and predicted yields and head weights for 
several Kansas wheat fields. 






Obs 

Pred 

Obs 

Pred 

County 

Field 

Variety 

Year 

Yield 

Yield 

Head 

Head 






2 

(gra/m 

) 


Riley 

ET Lab 1 

Centurk 

1975 

280 

431 

n.a. 

582 

Riley 

ET Lab 2 

Sage 

1975 

206 

439' 

n.a. 

593 

Riley 

ET Lab 3 

Tamu 

1975 

227 

426 

n.a. 

575 

Riley 

ET Lab 4 

Trison 

1975 

298 

436 

n.a. 

589 

Riley 

ET Lab 5 

Arthur 

1975 

297 

424 

n.a. 

573 

Riley 

Lysimeter 

Cloud 

1975' 

268 

436 

n.,a. 

589 

Riley 

Lysimeter Low 

Plainsman 

1976 

195 

3.76 

660 

511 

Riley 

ET Lab Low 

Centurk 

1976 

306 

343 

700 

463 

Riley 

ET Lab Med 

Centurk 

1976 

339 

436 

680 

590 

Riley 

ET Lab High 

Centurk 

1976 

289 

345 

680 

466 

Riley 

Cedarburg 

n.a. 

1976 

183 

254 

470 

343 

Riley 

Puett 

n.a. 

1976 

174 

292 

480 

394 

Riley 

T 1 

Triumph 

1977 

165 

153 

210 

207 

Riley 

T 2 

Triumph 

1977 

159 

209 

220 

282 

Riley 

T 3 

Triumph 

1977 ■ 

183 

174 

200 

235 

Riley 

T 4 

Triumph 

1977 

110 

160 

190 

217 

Riley 

E 1 

Eagle 

1977 

208 

210 

260 

280 

Riley 

E 2 

Eagle 

1977 

183 

252 

240 

341 

Riley 

E 3 

Eagle 

1977 

183 

200 

210 

270 

Riley 

E 4 

Eagle 

1977 

153 

161 

190 

217 

Riley 

P 1 

Plainsman 

1977 

128 

160' 

200 

216 

Riley 

P 2 

Plainsman 

1977 

147 

142 

220 

192 

Riley 

P 3 

Plainsman 

1977 

183 

143 

250 

193 

Riley 

P 4 

Plainsman 

1977 

165 

145 

220 

196 

Riley 

Cedarburg 

n.a. 

1977 

242 

372 

400 

*502 

Riley 

Puett 

n.a. 

1977 

164 • 

269 

250 

364 

Ellsworth 

E, Zeman 

n.a. 

1976 

176 

190 

n « d. o 

257 

Ellsworth 

A. Zeman 

n.a. 

1976 

302 

343 

n.a. 

463 

Ellsworth 

Bartunek 

n.a. 

1976 

335 

234 

n.a. 

316 

Ellsworth 

E. Zeman 

n.a. 

1977 

293 

290 

n.a. 

392 

Ellsworth 

A. Zeman 

n.a , 

1977 

116 

278 

n.a. 

375 

Ellsworth 

Wm, Bartunek 

n.a. 

1977 

331 

363 

n.a. 

490 

Colby 

19 

n.a. 

1976 

184 

261 

n.a. 

353 

Colby 

20 

n.a. 

1976 

268 

394 

n.a. 

532 

Colby 

29 

n.a. 

1976 

244 

375 

n.a. 

506 

Colby 

17 

n.a. 

1977 

159 

348 

n.a. 

470 

Colby 

19 

n.a. 

1977 

232 

306 

n.a. 

413 

Colby 

20 

n.a. 

1977 

226 

347 

n.a. 

468 

Finney 

18 

n.a. 

1977 

206* 

195 

419 

263 

Finney 

46 

n.a. 

1977 

167* 

193 

225 

261 

Finney 

50E 

n.a. 

1977 

153* 

335 

440 

453 

Finney 

129 

n.a . 

1977 

185* 

278 

481 

376 

Finney 

186 

n.a. 

1977 

240* 

199 

360 

269 


*SCS estimates 
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F±gi ?* Observed 'versus- pr-eddcted yiel’d for several' wheat crops-;. 
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reasonable estimates of dry-matter accumulation- during the winter. 

Equations- were developed- to estimate wheat yield from growth model 
input' and output data,. Reasons- for the poor yield prediction of the 
equations were discussed. 
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TRANSPIRATION-BASED WINTER WHEAT YIELD' MODELING IN DIVERSE ENVIRONMENTS 

ABSTRACT 

A modification of the longrstanding. evapotranspiration (ET) /yield 
relationship is proposed for winter wheat yield modeling. The model is 
based on output data from a comprehensive ET" model that has been developed 
and tested for several locations. Data included two years of yield and 
ET information- from sites across the state of Kansas.. The yield model 
perfo 3 nned favorably in predicting yields across Kansas for the 1977 crop 
year., for a--s±te- -in-Texas-,, and under conditions of induced diverse condi- 
tions by chlorocholine chloride (CCC) . Advantages- of this modeling approach 
include limited input data, simplicity of quantifying terms, and proven 
capability. Disadvantages include lessened accuracy at high, soil moisture 
levels, requires estimates of leaf area index and plant ontogeny. The 
model should prove- useful' where water' remains the- major limiting factor in 
wheat production. 
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TRANSPIRATION-BASED WINTER WHEAT YIELD* MODELING IN DIVERSE ENVIRONMENTS 

INTRODUCTION 

The ability to ascertain the yield of a crop by computer simulation 
(e.g. systems analysis) techniques has intrigued; agronomists, economists 
and agribusinessmen alike for many years. The- scientist is often fascin- 
ated by the end results of systems techniques, a better conceptual view 
of the systems, or mechanisms by which a plant determines its yield. The 
economist or agribusinessman is interested because crop prices (and asso- 
ciated market and futures profits) are directly affected by official- dec — 
larations of predicted world yields. Computer simulation,, coupled with 
remotely- sensed satellite imagery, offers a valuable tool to predict 
world yields- with- limited environmental information.. The. current world 
food and fiber situation with its Malthusian shortages, climatic stresses, 
etc. , give additional impetus to- the formulation of operational yield 
simulation systems. 

One method for predicting relative yields that has been given varied 
response for many years utilizes evapo transpiration/ yield (or more cor- 
rectly, transpiration/yield,, as- we shall show later) relationships. The 
direct linear relationship between yield and evapo transpiration (ET) was 
first comprehensively documented by Briggs and Shantz (1913). Since that 
time, the literature has -been replete with both supportive and skeptical 
reviews. A few of those that are directly related to this report are the 
positive aspects shown by DeWit (1958) , Arkley (1963) , Stewart (197'2) , and 
Hanks (1974) . Rawitz (1970) showed data that indicated the relationship 
between ET or T (transpiration) and yield is not completely linear or 
correlated, especially at high water levels. There are others who question 
the relationship vigorously, but this is no.t well documented in the litera- 


ture. 
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We currently have a computerized ET model at the. Evapotranspiration 
Laboratory that is operational for a number of sites across Kansas and 
the Great Plains. The model is designed so as to accept satellite sensed, 
crop data (leaf area index, LAI)-, LAI measurements- give the model a dy- 
namic response to crop growth and allow the splitting of ET into its in- 
dividual components of transpiration (T) and soil evaporation (E) . 

We wished to kno.w whether winter wheat yields could be predicted with 
reasonable accuracy, utilizing the comprehensive output data of this ET 
model. A yield, estimate based upon physical quantities- and simulated 
crop responses was preferred over a purely statistical approach that is 
not truely dynamic in its response to changing crop conditions-. The ques- 
tion to be answered was simple: we have the ET data — can it be used to 

predict yields over wide areas and diverse environments? 

THEORY 

Evapotranspiration and- Yield — Dry Matter 

The rationale- of ET based yield models of the past, has been that the 
processes of biomass synthesis and translocation are closely related to the 
transpiration pathways within the plant (Boyer and McPherson, 1975), thus, 
we exp.ect that translocation cannot proceed, without transpirative flow and 
that the metabolic rates would follow the transpirative response to water 
stress. Further, we anticipate that the critical processes of photosynthesis 
and CO^ gas exchange will be directly correlated, to rates of water vapor loss, 
at the stomates,, as D'eWit (1958) suggested. The correlations are indeed- 
complex and not fully linear; however, the large portion of reported data 
would be supportive that yield is highly correlated to plant-soil-water 
relations. To the extent that water is limiting, T/Yield relationships 
must be employed to adequately predict yield. In Kansas-, and throughout 
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4 




much of the world, water relations are the major limiting factors in 
crop production. 

A mathematically-sound expression of the transpiration/yield rela- 
tionship was given by Hanks (1974) as: 

Y/Y = T/T [1] 

P P 

where Y is actual dry matter yield, Y is the dry matter yield when water 

P 

is not limiting, T is the actual transpiration sum, and is the poten- 
tial transpiration sum if soil water was not limiting-. Note that any lim- 
iting factors other than water will be acco.unted for in the Y^ term that 
must be determined for the given site conditions and year. This iS' a 
simple expression of, the direct limiting influence of water stress- on 
yield. 

This equation has applicability only under conditions of differential 

water stress and- is not intended to model any or all of the many variables 

that can be included in the Y term. Hill, et al. (1974) showed the value 

P 

of this approach for irrigation management problems-. 

Most investigators, previous to Hanks (1974), correlated yield with 
ET rather that T because in actual field conditions, ET. cannot be measured 
separately as E and T. Because soil evaporation is seldom constant or 
negligible in grain crops, the combining of E with T must surely introduce 
an error term into the yield relationship-. With energy balance, the 
evaporation from the soil surface and from the plant can be evaluated sepa- 
rately (Ritchie, 1972)'; thus, equation- (1) can be evaluated with reasonable 
accuracy. 

Evapo transpiration and Yield — Grain 

The processes that govern grain forma-tion within the plant are much 
more complex than those producing dry matter (biomass). For example, a 
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wheat crop may have a lush vegetative appearance, yet have little grain 
at harvest. It is- agreed that in determinate crops, such as wheat, certain 
stages of growth are more crucial in transpiration/yield relations than 
others (e.g. flowering and grain filling). This has been documented in 
excellent reviews such as that of Boyef and- MacEherson (1975). Hanks 
(1974) revised an equation of Jensen (1968) to include relative transpir- 
ation, rather than relative evapotranspiration: 

n , . 

Y/y- = TL (T/T ) ^ [2] 

P i=l P 

where Y is actual grain yield, Y^ is grain yield if wate r is. n ot- lim iting, 
i is a phenological stage of growth, T is actual transpiration for that 
stage of growth, is potential transpiration for that stage if water is 
not limiting., and A is an exponential weighting factor. 

Equation [2J has been shown accura.te for corn (Hanks, 1974; Hill et 
al. , 1974) and spring wheat (Rasmussen, 1976) under conditions- of differ- 
ential water stress; therefore, it is an extremely useful approach in 
irrigation management studies. However, Y^ must be determined for each 
variety and locale on a yearly basis especially if the length of the pheno- 
logical stages vary. 

Evapotranspiration Model 

Kanemasu, -ef al. (1976) described a versatile ET model that has been 
developed in computerized, form. The model uses, inputs of daily solar 
radiation, LAI, maximum temperature, minimum temperature, and wa.ter re- 
ceived (precipitation and irrigation) , to compute daily values of soil 
evaporation (E) , transpiration (T) , potential E and T, and a detailed 
ledger of the water balance within the soil profile. Season inputs are 
soil profile characteristics- such as water holding capacity and surface 
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dirying- characteristics.. This model has been shown to be- effective for 
soybean and sorghum (Kanemasu,, et al., 19.76) wheat (Kanemasu, et al., 1977), 
and corn (Rosenthal, 1977). 

This model is initiated with a layered soil profile of known water 
content and water holding/loss- characteristics-. Each day the- energy avail- 
able for water loss is evaluated' with the. modified Priestley-T-aylor equa- 
tion as described by Tanner and Jury (197'6).. Soil evaporation and crop- 
transpiration are evaluated- separately.. Surface soil evaporation is 
allowed to proceed at its potential rate during the day following an addi- 
tion of water (rain or irrigation) . After' the soil has lost a certain 
amount of water, E decreases with the square root of time until, the soil 
reaches an air dry lower limit or more water is- added. Transpiration 
proceeds at a potential rate until 65% of the- soil's available water is 
depleted, and, then decreases linearly to- zero at the permanent wilting 
percentage (PWP).. The energy is partitioned between the evaporation- and 
transpiration processes according to the LAI. 
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YIELD' MODEL DEVELOPMENT 

We desired a model that would respond to the differing, conditions of 
the growing season throughout the Great Plains. Differential water stress 
is the most likely limiting factor to wheat yield in this wide area. How- 
ever we realized that water is not often limiting in the eastern belt of 
this region where annual precipitation approaches. 100 cm. A small data 
set from throughout Kansas for two years (see Table 1) was assembled for 
model formulation. 

The growing season was separated into three phenological periods; 
emergence to jointing (approximately October to March in Kansas), jointing 
to heading (approx. March to May), and, heading to soft dough (approx. May 
to June) . This, separation allowed the period of rapid vegetative growth 
(jointing to heading) and rapid grain filling (heading to soft dough) to 
be evaluated independently. Phenological devlopment was modeled using 
a modified fom of the Robertson Biometeological Time Scale (Eeyerherm 
and Paulsen, 1976). 

The output of the ET model was programmed to show the following quan- 
tities evaluated for each phenological stage and for the entire season 
(see Table 2) . These summational relationships will be- referred to by 
number as model development is discussed. 

Hanks (1974) and Hill et al. (I974-) used relationship #4 in Table 2, 
in their multiplicative grain model for corn. However, the X values, were 
determined by empirical (trial and error and experience) means. We desired 
a more quantitative approach. Regress ion- type models were readily avail- 
ble as statistical packages- for computers and programmable calculators; 
thus, initially a simple additive multiple linear regression model was 
tested. The square of the Pearson product-moment correlation (R ) for 
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Table 1. Summary of components- of calibration data set. 


Calibration data set 
1974-1975 


Ashland,' Riley County, Kansas 

5 varieties on- intensive E..T. Lab site 
(Centurk, Arthur, Sage, Tamu, & Trison; 
with yields- averaged to one value per variety) 

3 cooperating farm- sites (one ave. yield value per site) 
Ellsworth County,, Kansas 

3 cooperating farm sites (one- ave. yield- value per- site) 
Finney County, Kansas 

3 cooperating farm sites (one ave,. yield' value per site) 
TOTAL: 14 observations 


1975-1976 

Ashland, Riley County, Khnsas 

3 planting densities- of Plainsman V variety 
(one ave. yield value per density) 

3 cooperating farm sites (one ave. yield value per site) 
Thomas County, Kansas 

3 cooperating farm, sites (one ave. yield value' per site) 
Ellsworth County, Kansas 

3 cooperating farm sites (one ave. yield value per site) 


TOTAL: 12' observations 
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Table- 2 . Normalized actual/, potential .wa'ter- loss fractions used:- 

1 ,. s:(.T/Tp) 

2. S(T/ET ) 

P 

3. r(_ET-/ETp)- 

4- , (ZT)/(2Tp) 

5- CET):/(ZETp), 

6. aET)-/‘(SET^p) 

Where-: Tj = Actual- transpiration; (mm/ day) 

T^=- Potential- transpiration (mm/day) 

ET= Actual .-evapo transpiration,- (mm/ day) 

ET- = Poitenti'al'.' evapotranspiration- (imn/'day)’' 

P 

All values- are daily values. All, summations are for 
a given phenological- stage. 
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observed versus predicted data points was used to evaluate the suitability 
of each model formulation. 

The simple additive regression, model of the form-: 

Yield = Constant^ + C^(X1/X1') + C^CXY/XZ’) + C^CXS/XS') [3] 

where Xi/Xi is one of' the six ET relationships- in table 2, gave first order 

2 

results with formulations #1 and //2’ (R ^ .4'). However, this formulation, 
though statistically sound, did not seem reasonable in the physical sense. 
It is not reasonable to assume that with complete deficits in stages two 
and three, a grain yield could still be achieved. This model was- abandoned 
in favor or the more reasonable Hanks/Jensen form. 

By using a logarithmic transform: 

Log(Y)= (ConstantQ)+ X^Log(Xl/Xl’ )+ X2l'OgCX2/X2' )+ X^Log(X3/X3' ) [4} 

with a simple multiple regression routine, the gamma values for each stage 

of the Hanks/ Jensen equa-cion can be determined statistically with Y^ for 

that set of data statistically lumped into the leading constant C^. 

Utilizing this technique,, each of the ET relationships in Table 2 were 

evaluated in the Jensen/Hanks grain equation. ET relationship #'4, which 

is the sane- as Hanks (1974-) and Hill et al. (1974) used for corn and 

2 

Rasmussen (1976) used for spring wheat, gave poor results (R ^ 0.2). 

Initial inspection- of the data revealed that for most of the data, water 

stress did not occur in the- test sites (most were in Eastern Kansas) . 

Therefore,, yield differences could not be accounted for. 

Further modelling using ET relationships #1 and #2 gave much better 
2 

results (R - 0.4 to 0.6). Nvimber 2 gave the best overall results through- 
out the entire data set with the following equation showing the highest 

2 

statistical significance (R “ 0.6): 
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Yield" (■■bu/acre) =• 

2.856: (E(T/-EET)-)^’^°'^ (X('T/EET,).) ^ [5] 

Where a=growt-h stage fr,om emergence, to, jointing- . 

B=growth- stage from' jointing- to. heading 

c=growt-h. stage from heading, to- soft-dough 

The observed: versus:--predicfed yield data, .are piotted' in -Figure 1. 

The- l'9'7-4— 75' data co.vers-. a- broad range of yields and the- correlation is., high 

(-.80,). The-, 19-7 5--7.6i data- set, showe.d extreme variations- within samples.. Corr 

relation with .this data and modeled response: is low.' ('.56).., The’ pooled 

2 

correlation is- 0..7.3 or- an-R- of, 0.54-i 
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MODEL EVALUATION 

Additional yield data was- obtained for* the crop' year 1976-1977 for 
sites across Kansas. In addition, a data set was obtained, from Bushland, 
Texas in which the required LAI and soil moisture parameters were avail- 
able- The Kansas data included data from an intensive- experiment in. 
which CCC (a plant growth regulant) was applied to three varieties to 
induce variate -leaf- area and water use conditions-. Considering the wide 
variation of- annual rainfall across Kansas and. the variate conditions 
induced, with CCC, a diverse set of test conditions was imposed. 

A summary of these yield- data as compared with modeled yields is 
given in Table 3, and plotted, in, Fig. 2. The pooled Bear son-moment 
coefficient is 0.81. 

Additional attempts were made- to modify the model by inclusion of tem- 
perature- and potential yield, terms-; however, no significant increase- in 
Pears on-moment was ob.tained when using the entire data- set. 
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T-able 3. Simmia-ry of observed versus predicted yields for all yearsr.. 


KS 1975 

OBS PRE OBS/PRE 

KS 1976 

OBS PRE OBS/PRE 

RS 1977 

OBS PRE OBS/PRE. 

BUSHLAIR), TX' 
OBS PRE OBS'./P-RE 

- 

45- 

- 

30 

31 

.968 

34 

32 

1.0.63. 

55 

46 

1.196 

29 

38- 

.763 

44 

40 

1.100 

30 

29 

1.-034 

40- 

40 

1.000' 

20 

35 

.571 

40 

36 

1.111 

30 

30 

1.000 

50 

43 

1.163 

32 

40 

.800 

50' 

3,9 

1 .. 282 

25 

25 

1.000 

45 

44 

1.023 

49 

47 

1.043 

53 

52 

1..058 

21 

30 

..700 

40 

42- 

.952 

49 

45 ■ 

1.089 

47 

45 

1.044 

24 

26 

.923 

40 

41 

.976 

37' 

• 44 

.841 

29 

37 

.784. 

30. 

26 

1.154 

45 

41 

1.098 

34 

34 

1.000 

- 

40 

- 

27 

26 

1.038 

20 

22 

.909 

46 

46 

1.000 

- 

47 

- 

27 

25 

1.080 

13 

21 

.6-19 

24 

25 

.960 

55 

36- 

1.528 

26 

29 

.897 

35 

38 

.921 

20 

30 

.667 

29 

31 

.935 

30 

26' 

1.154 

29^- 

30 

.967 

39 

34- 

1.147 

49 

40 

1.225 

18 

26 

.692 

42 

44 

.955 

25 

19 

1.316 

23 

39 

.590 

19 

31 

.613 




17 

16- 

1.063 




48 

33 

1.450. 




24 

17 

1.412 




54 

44 

L. 230 










27 

34 

.794 










40 

38 

1.053 










26 

33. 

.-788 










38 

28 

1.357 










37 

36 

1.028 





Pooled Corr. Coeff. 


.81 
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DISCUSSION 

Emphasis- was made,, at the beginning of this paper on the need for an 
operational yield modeling system for winter'wheat, in- Kansas and other 
areas of the world. We believe that the model discussed here will provide 
first order results under conditions where water is the major limiting 
growth factor. No attempt is made in this- model to consider the effects 
of fertilizer management problems., tillage practices, etc. However, con- 
sidering the fact that yields can be modeled with reasonable accuracy across 
the entire: state of Kansas and in Texas,, we would again propose that water 
limits yield to a considerable extent in these areas.. It is evident from 
this study that has been performed that other factors than, those controlled 
by model inputs (temperature, radiation,, rainfall, soil water character- 
istics, etc.) need to be accounted for in future modeling endeavors. 

Considerable stress was. placed upon the fact that relationships 

should be the major indicator of water stress'. The model formulation- 
that gave the greatest degree of accuracy was the T/ET.^ form. There are 
several reasons for the statistical suitability of this term. Water was 
not limiting to any extent in all data sets from Ellsworth and Riley 
counties. However, we were modeling winter wheat which has a long period 
of dormancy in the winter. According to the model f/Ip would approach 
one for this entire period (even though T^ is near zero-)-. However, if 
T/ETp is used,, this number reaminas low — being an indicator of little 
growth or photosynthesis. 

The summation form is used (#1 and :/2) instead of the fractional form 
(//4 and //5) , due to increased statistical credibility because a summation 
form allows a given periods contribution to yield to be larger as time 
within the growth period is lengthened. Thus, if cool nights produce a 
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longer, grain filling- period in Western- Kansas.-, (as- they often* do)', yield' 
is allowed to increase- The-, fallibility of this* approach lies, in its 
■physical' meaning' — that- yield* could* approach, inf inity if grain; filling, 
were prolonged, to infinity i Evidence suggests; that, there is a genetically 
determined upper limit of., pheno logical periods*. ' 

In. reality then, the* present model [.5] integrates, a. number" of - compo- 
nents of yield -into one* equation: . rate of growth (as indicated by T/ET ), 

P‘ 

length of pheno logical period (culmulative photosynthesis) , and water 
stress. Further refinements are currently being investigated- The- indi- 
vidual .components are. best defined* if isolated. — not lumped,, as in. the* 
present, model; however, the model, can predict .yields well, within the* 
range of", field variab-ility and, should provide reasonab'le yield estimates- 


where water- is- limiting-. 
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AP'EENDIX'. A 

ALGORITHM. FOR DAILY GROWTH -MODEL 
FOR. WINTER 'WHEAT 
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KSU HCDEt 6: SOIL gViPCRATICh HCCEL. 

INPUT CAROS. 

1. TITLE CARO. APPEARS AT THE TOP CF EACH 
OUTPUT PAGE. 80 COLUMNS OF INFORMATION. 

2. CROP DATA. OATE SPECIFIED AS OCT 4. 1974 IN FIELD OF 32, 
LOCATICA, lA FIELD OF 32. 12 SEPARATE CAROS). (4A4), 

3. VARIABLE FORMAT. SPECIFIES FORMAT CF INPUT 
FCR HCCEL. INPUT CARD *7. 

4.. PARAMETERS FOR RUN. FIELDS OF SFIO. 0.6F5. 0, 

5. THE HCKIH DAY AND YEAR CF START CF DATA.. 

SPECIFIED AS 02/14/75. {1 2 .IX ,{ 2, IX, 12) 

6. THE LAST HOAt). OF DATA SPECIFIED, 12. 

7. CFLAG FCR SCRGHUH=l ; S0YB£AA=2; HHEAT=>3j CORNEA .. (12) - 

8. THETA VALUES FCR EACH LAYER. THETA HAX VALUES FOR EACH 
LAYER. THETA MIN VALUES FCR EACH LAYER, (1SF5.0) 

9. CHECK FOR 8HTS SUBROUTINE. 01 IF NO, 02 IF YES.. 

CHECK FCR IRRIGATION. 00 XF NC., » CF DAYS ♦ 1 IF YES. 

10. STAGEA CAYS FROM EMERGENCE. STAGES DAYS FROM EMERGENCE- ' 

e.G. 34 59 FORMAT (13., IX, 13) . 

11. COEFFICIENTS USED IN BIO HETERLOGICAL TIME SCALE 
ISU8RQLTIAE tCLGKER)). 6 CAROS, FOR VHEAT ONLY. 

12- THE MCNTH ANO DAY OF PLANTING. FOR WHEAT ONLY. (12, IX, 
12 ) 

13. DATACARDS: OA I LY SR, M AXT ,H IN T, LA I, RA IN ,0L 

SHOULD CCRRESPCND TO VARIA3LE FORMAT #3 ABOVE, 

OL FOR WHEAT CNLY, 


DECLARATION OF VARIABLES 

REAL MXH2C, HAXT, mint , LAI .THEVAL tSJ , T VAL ( 5 ) , ZVAL (5 ) /50- , 250 

• 2*300., 900. /.CVAL ( 5 1 ,TFE MAXI S ) , THE M IM 5) , KVAL2T 5) / . 05 , . 25 , .25, 

• .25,,2/-.KVAH(5)/. 1, .5 ,.4,2*0./,TH£TT.TAU, 

•RAIMOU.RAINEH.TACC ,HLLT,TDAY,TX,TN,0L, 

• MFLDCF ,KS ,T , £T . FUOCP. T2. CV, 

• PTACCl 

REAL L3I7)/'PT -•.•EM -'.'JT -< ,»H0 - • , • SO -•,'RP->*/', 

*O0(7)/'EH* jr ' , •6T’ FD', 'SO* , 'RP' , ‘ »/,T3 
REAL KA7( 61/1-0.2- a. 2. 7, 3, 0 .4. 0 ,5. 0/ .0 1 FF .BTST L ( 7, 13 )/9 1*0.0/, 
•T0TL{6)/6*0.0/.BT0TLI 13 )/13+0,0/ 
integer SYR(6>. SMC(6).SDAY(6 1 .SCLK ( 6 ). CLOCK, 8TFLG..IC0UNT.RFLG, YEAR 
REAL LVAL(5)/50.,2 50..2»300. ,600./,MCNrL( 12.14) 
real INTL.GRGSPH.CSESP.NRSSP .NCE.OMTR, DH.CDH, TRcSP ,DLAI 
INTEGER CAY.YR.SCYSCR,CAL(12 ) /3l .2 8 , 31 . 30 .31 ,30 .31 , 31 , 30 ,3 1, 30 ,31 
*/.TlTLE(20) .FLAG.CCI 5) /5*0/,0AYT,PLANTl8) ,LOC( 10) . 

•K.POAY .PMC.BMT.efLAC.CTR.J ,DT .N.OHS 


S 
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REAL PLAJ,LAlHAX,TEO,ETEC ' 

JNTECER FL-AGL. START 
IMTEGER SGSA.SGSB 

integer IFCOOI.IOAYtSOl ,IYR(30) .IRRCHK.I I 

REAL. IRR130) tFALLBT, FNTE><P- 

INTEGER IFLAG 

PEAL IBRIG 

DIHENSICA CQEFlfi.SJ 

REAL*a FRHTMOl 

initialize 


II=«1 

RAINCL=>0.0 
EST=0.0 
£T=0.0 
CCUNT=2.0 
CAYT.=0 
TACC=0.0- 
K=l.. 

KS=l.O 
aHT=0.0, 

T2=0.0 
TAU=0-0 
Cy=3'. O’ 

CT=60. 

FT ACC 1=0,0 
5TFlG=0 
RFL6=0 
CLCCK=0 
ICOUKr=0 
FriTEHP=15.0’ 

CO’ 20A J=l,6 
SYR (J 1=0 
SHC I-J 1=0 
SDAYI J)=0 
20^ SCLX(JJ = 0 
START=0 
INTL=0,0 
GROSPF=Q.O 
CR£SP=0’.0 
NR£SP=C.O 
TRESP=0.0 
AGE=0.0’ 
CKTR=0-0 
CH=0.0 
COF=0.0 
FLAI=0.0 
FLAGL=0 
LAIHAX=1.0 
FTP-m=15.0 
CO' 202 J*l,l2 
CO 203 I = I.IA- 
203 PCNTLtJ.I 1=0.0. 
202 CONTINUE 

CO 210 J=1.7 
CO 211 I =1.5 
211 STSTHJ, I)=Q.Q 
210 CCNTINUE 
C 


ORIGINAL PAGl Ife 
OF POOR QUALi^re- 
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C READ tN INPUT 

C CUT.PUT HEAD tNGS. 

C 

READ- t 5.-100. EN0=S5?J- TITLE' 

BEAOt-S.lOOJ PLAhT 
REAOIS.lOOl LCC 
REAOIS.ioriFRHT 

RtAO' (5.1021 PXH2C.TV.TIN.T5 .U.ALP HA.FALLfiT.CNSJ.X5.FL OCP .HOLT 

TFETT»TIN 

KFLOCF=FLCCP*1'500> 

WRITE ('6.1Q3)TITUe 

WRITE 16, 10A-)HXH20,TV, TIN.TS, U. ALPHA .F-ALLRT, CNST, X5 ,'HFLDCP »HULT 
WRITE! 6 „r05) PLANT 
WRITE (-Oil 061. LCC 
TI.TIN 

READ (S'.nni HO, DAY, JR 

YEAR-YR- 

KOAY-CAY. 

read; (s; 1C7T KHO^ 

IF(KKC.LT.K0]XH0=KH0 + 12’ 

KHO=HC ' 

HMH=HC. 

PEA0(5,107)S0YSQR- 
ALPHAV=1',A4 • 

GO TO. {301,3a2.303,3C41<,.SQySOR. 

302' AUPHAV=1.71 
GO TO’ 30 r 
303 AtPKAV=1.56 
CV=l.3S 
0T = T'S5. 

GC TO 301 
304'-ALPHAV=li74" 

301 READt-S'.lCaT THEVAL .THEFAX.THEHIN- 

00-200 r=i,s 

200 CVALl.n=(-THEHAXn )*.,3H-(THEHIN(n*.7,) 

READ (S..1071 BFLAG.IRRCHK 
IFISCYSCR.NE-.U GC TC 501 
REAOi 5.121) SGSA.SGSB 
501. IF (8FLAG-2); 305.306,305 
306 READl.S.lCS) -( (CCEF (T , J ):, J=l ,8}., I=V, 6 1'- 
REAO (S.-lOTl PPOiPO'AY 
305 IFURRCHKiEO.O) GC TC 359 

ReA0(5,l'20HIHC( I J , ICAVni , I YR ( I ) , I RR( 1 1 . I = 1>.1RRCHK) 

J=IRRCHK-1 
WRITE (6,. 1221 

HRITE(6.1231 (>IMG! I) . I DAY ( I )' ,I YR( M , IRRI IT, I J) 

359 FLAG=0- 

Ifi (T5.-GE.X51FLAG=l 

TOP OF-UCCP OF ANALYSIS'.- 

206. WRITE (6.103 1- TITCE. 
kRITE(6.llO) 

uf=cauffci- 

IF(HMC.EO,2,AND.HOO( (YR>1> .4},, EQ.O) LIM=2,9 
DO; 207 JJJ=KCAY.L1K' 

IFIHNC.EO. l.AND.JJ J.EO.l.) YEAR=,YEAR 

IFtSCYSCRlEO.lT GC TG 308- 
GO TO. ('307V3081 .BFLAG ' 




OT 
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308 Re’A0’(5,FR«T,ENC=332) -SS‘, HAXT.'.HINTtL AURAINiOL. 
GO -TOl 30<5- 

307 ■ READ ( 5 i.FRMT , EN0=332 ) SR , RAX^t H I NT", LA I:. RAIN. 

309 COveR'=LAI7CV> 

IFICOVER.GT.0.0) .0 AY,I=5C AiT.T+ 1. , 

IF(CDVER.GT.I.O) CCVER=‘1'.00 

IF JOAOrT-O.T) -333.333,310,' 

310 IF(CCVER-.LT-..4.1 cover*. 4; 

IF tehpurature; is in- cegrees- f;,. 

PLACE' CONVERSIONS f=t;T-32'. 0 1*5. /9i„ FOR' 
HIM AKO-MAX~hERE. 


333 •FAXT-=n'AXT-32-.Ql*5‘i/9.. 

Fi NT= ( KI‘AT-32. 01 *-5./-S. 

1FLAG=0- 

IF'MRRCHF'.EO'.O) GC TC'36'l 

IF-n IF.aUl ).£OCHHa.ANO'.,ICAY( III-.EQ'.JJJ) .ANOS-I-YR II I J-Ea.'.YEAR 1 
*GO TO; 363' 

• GO- TO' 361. 

363 RAIN=RAlN-t-lRR'{:in” 

IF.L'AG=1 ' 

IR8IG=.IRRIin*I0,. 

I r=.I5['+l 

361 RATN=RAI'N*lQ-.0' 

RAIN£V.*RA1N>RAINGL 
RA !N0L=RAIN' 

IF' IRAINEH.LT..6.0}' ,GC TO 3tl, 

FL'AG'aO 

EST=0.0 

CCUNT-=2.0 

31i THP*L3*HAXT.+HINT-T/4.0: 

■ SSO^OELT'AITNP) 

313: GO' TO- t3l2.,315),9FLAG 

315 IF IMHO'.EO,.pRO,ANO*.J JJ.EC.POAY l. SHT=l ' 

IF (BMT-U.312:,316,3 12 

FCR WHEAT- ONLY- 


316' T.X*MAXT ‘ 
TN=HINT 


CAUL.-BiO-HETERLCGICAU TIHEl SCALE SU6RCUT-INE: 

CALL- CLQKER' ICGEF.-TN.-TX , CL ..TCAY-, X-. KU LT, T ACC. CLOCK, HHOrJJJ .-YEAR , 

* BTFLG.OIFF ,SYR’,SHO.-SOAY-.SCLKVICGUNTSRFUG,.KAX'» 

♦ PTACCl) 


CALL CF SUBROUTINE. TO' CALCUliATE' POT EVAP.. 
312 CALL PCTEVA(LA'I-,SOYSOR.RN.SR,ALPHA',-SSDs£OVOA-Y.T1 ' 


TI=TI-£TtRAlN 

IFITI.GE.PFLOCPl'. T.1 = )'FLCCP' 

C; 

C IF POT. EVAP IS- . ZERO' ALL. EVAP ARE' ZERO' 
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c 

IF teC.GT.O.OI GO TO 317 
60=0.0 
es=o-o 
T >0.0 
T2=0.0 
A=0.0 
GO TO 325 


CALL TRANSPIRATION SUBROUTINe 

317 CALL TRANS(TI,HXH2ciTV.KS,LA'I tT , ALP FA, RN, SSO, COVER . TAU , ALPHAV, 
♦ DRY,SCYSaR.T2,TGCQ‘,TACC, FALLRTJ 


CALL SCIL EVAPORATICN SU0RCUT.INE. 

CALL EVAP(FLAG,LAI.,RK,EST,eS,U,CNST,COUNT,SSO, DRY, COVER, OAYT.TAU, 
♦SOYSQR) 


CALCULATICN OF ^ EVAPORATION 
CO TO 1315,320,321,319) .SOYSOa 
FOR SOYBEAN 

320 IF lHAXT-31.) 322.318,318 
FOR HFEAT 

321 IF {HAXT-27.) 322,323.323 
FCR SCRGKUH OR- CCRN 

3A9 IF5MAXT.CE. 33,0) GO TO 323 

322 A=0.0 
GO, TO 32A 

318 A=.25*T 
GO TO 32A 

323 A=.1*T 

324 lFIHAXT,.GE.-3.0) GO' TO 325 
ES=0.0 
A=0.0 

325 IF5KS.LT.1.0) A=0.0 
T2=T2+A 

TAVAI,L>TFETT-ITVtl500) 
KS=TAVAU/(FALLRT*MXF2C) 

IF (KS.GT.U !tS=l 
IF IKS. LT. 0.0) KS=0.0‘ 

T3=T2»KS 
ET=T3*ES 
C 

CALL 0 ISTR (-T3 , T VAL , KVA L-l , K VA L2 , L A I ) 

c ..... 

RUNOFF=0.0 

0RAIN=0.0“ 

CALL CAYl (THEVAL,ZVAL , DRAIN, CC,THEPAX ), 
IF (RAIN.EO.O.)GO TO 401 


ORIGINAL PAGE IS 
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IF('IFL;SG.£0.il} Ca TO 3'66 
GO TC 367 

3<6 BAI^l = PA^^-It^RIG 

CALL CATC( TrieVAL.RAINrZVAL.RUNCFFtCC.THeMAXJ 
CALL CAYZtTFEVAL. IRk IG.ZVAL, RUNOFF .CC.THERAX) 

rain=ka:s+irrig 

GO TO AOl 

367 CALL CAYOlTHeVAL.RAIN-.-ZVAL.RUNOFFi CCtTHeMAX) 


401 CALL RCIST ITH6VAL . £S. TVAL . ZVAL .THEM INO 


IF( EO. 68-0,0) GO TO 405 

T£Q=0.0 

ETEO=C.O 

GO TO 406 

405 T£C=T3VS0 
£TEQ=ET/EC 

406 THETT=0 

DO 201 Jal.S 

201 TKETTitTFEVALl JjALVAUJ) 1 + THET.T- 
GOTO L407,402,403» lOOO 1 ,SCYSOR 

402 CALL JIHILAI ,PLAI,.INTL .GRCSPH.-TRESP iNCEf DF.CDHtSRtSTART.KS ) 

GO TO 404 

403 CLAI =LA I + ,02'5. 

IFtOLAI.GT.4.62) DLA.I=4.42- 
INTL=5.21ASR*0LAI**.32S6 

TCP'S PBCGRAH 

CALL TCMtFLAGLrLAI ,PLA I-KiTACC-.GROSPH.SR', PAXT.ORESPiCDH, 

♦ DL.NR£SP.NC£.'OM,«INT,PTH'lN,LAlPAX,K-S) 

TReSP=CReSP*-NRESF 
GO TO 404 

407 CALL SCFGPH (CAYT-.HAXT.H IhT, PTH’IN ,SGS A, SGS8 . [NTL , SR , LAI ,GRQSPH, 
*DRESP .EL ,COR,NCE,KR£SP .OJ'.PLA-I .LAIMAX. FL A GL . PNI EHP', KS 1 

TRESP=CR£SPtMRESP' 

404’ IF(-T3.£Q-0,0J; GC TC 335 
DMTR=CP/T3 
GO TO 336 

335 CNTR = C'.0 

336 PTHU = PINT 
FLAI=LAI 


PRINTOUT CAILY VALUES' 

WRITEI6, 111)»1N0.JJJ,SR.MAXT,‘'INT,RAI6,LAI ,T3,£e',ET, INTL , 
*GRCSPH,CRSSP‘.NRESP.T'R£SP tNC’E ,CHTR, OPvCOHt TACC,-ThETT,KS 

PCNTL IPHQ, n=HCNTL(HPO. I J+SR 
HCNTLIPH0.2) =HCNTL(PPC.2 )»RA'1N 
KQNTLI FMG..3)=MCNTLIKK0 ,3}*-T3 
NCNTL IFFOt4)=HCML(P.''0.4J*£0 
P0NTL(FMC.5)=HCNTLLMP0.51 + ET 
MCNTUFHC,6)=HC6TL(HF0,6 )»INT-L 
PCNTL(FPC-7)=MCNTL(PFQ .7J+GRCSPH' 

HCNTL (FPC.8)=KCATLIPFO,a-)tORESP 
HOST L( PMC. 9) =PCML(«.‘'Q.9 )cNR£SP 
PCNTLtPMG, 10.1 = KGNTL( PMC . 10) + TRESP 
HCNTUFPO.M ) = PCNTUPHC, U J+NCE 
HONTLiHHO, 12 )=HC6 T L I FHC , 12 )F DHTR 



Ma^TU)'«0,^3l = t'C^TL(^'HC. 13J»CM 
HCNTUI'I'C, l'»)=Ci5M 
8HS=K 

IFIRFLG.ee. n. BFS=BMSvl 
8TSTLISHS, l)=BTSTL(ByS ,n ^GRCSPH 
BTSTL (E“S,2) = eTSTL ('BFS . 2 1 + TReSP 
BTSTL(3FS,3»=8TSTL{BFS.3 ) + NCE 
eTSTL(SHS,V)=8TSTLt3FS.'tJ + 0KTR 
BTSTLteJ'S,5)=CCF' 

8TSTL ( BMS ,6) =B TSTL CBFS .6 1 >T3 
BTSTL(e''S,7l=8TSTLIBFS.71+6T 
6TSTL (SHS.B) =3TSTUB)'S .8 1+160 
BTSTL t SHS,-;)=STSTL(BI<S .Bl+STEC 
BTSTL<8.“S.101 = 8TSTLI BUS. 101+KS 
8T STL(3FS.ll )=BTSTL.( BHS.111+T2 
8TSTL{BHS,-l2)=BTSTL(eHS,.12).+ SR 
BTSTLf BHS ,l;3)=8rSTLJ BMS.13 1 + INTL 
207 CONTINLS 

332 V.RITe(6.U21 1KCML{FH0 + J1,J=1.13) 

KHG=HHC+l 

IFtHHC.GT.12-J FMC=HHC-12 
KDAY=1 

IF tHKF.EC.KHOJ GC TC ?99 
FMF=KHH+t 
GO TG 206 

■999 kRITEte.lOSJ TITLE 
hRITEt6,ll3). 

HHO=HC 

DO 212 xy»>=KC.KHC 
VRITEti.mi HHC,«CN1L (HHQ,2) .HQN’TUHHC.Sl.HCNTLIHHC.Tl , 
=fHCKTLtHHC,101 .HOHTLtFHC, 1 1) , HGNTU HHO, 12) .KONTU HHO , 14) 
TOTLfl l=TCTHl)+FCRTLtHHC.2) 

' TOTU 2)=TCTU2J+HCNTL(HMC.51 
T0TH3 l=TCTLl3)+t'CNTL(HMC,7) 

T0TLt41=TCTL(4) +JCNTU HHC.IO) 

TOTL ( 5 )=TOTL (5) +HONTU HHO, 11.) 

T0TL(6)=TCTL(6) t!'CKTLtHMC.12) 

HMO=HHO+1 

IF1HHC.GT.12.) HHO=l 

212 CONTINUE 

WRITEt6.U51 {TCTL(J).J=1,6J 
IFtBFL3G-2J 1000..998 , 1 OQO 
998 WRIT£.{6.116) 

00 213 J=l..eMS 

Vft ITS (6,117) LSIJJ.UBU) , 1 8TST LI J, H ) . H=l , 13) 

CO 214 H=sl,l3 

214 BTCTL IH) = ET0TLIH)+8TSTL( J.H) 

213 CCNTINOE 

VRITEtB.llB) {STCTL(H),H = 1,4) , ( B fOTL I H) ,H=6 ,13 ) 

WRITE (6. 1 19) UAYIJ j.SCLKtJ,) .SMOI J) ,SOAY(J) ,SYRt J) iJ.= l.K) 

100 FORMAT t20A4) 

101 format- UCA8) 

102 FORMAT (5F10.0.6F5 .0 ) 

103 FORMAT ( 'I',20X.20A4.///) 

104 FORMAT I'Qi . 19X. 'HA^iiyuK AVAILABLE HATER tHH],., 

*• F10.4.//.20X. ' ThETi SUB V (15 EAR) 

* F10.4.//.23X. 'THETA INITIAL IN 5 FT. PRCFILE (HH) 

* F10.4.//,20X,'TFETA SUE 5 CM. LAYER.... 
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+ SRESP'rOH'.PLAI,l:AmA'X ,KS' 

RfcAL CTEJ'P.N'TSPP.Lf-.-TSTRtPNTEHP’.LMUI 

IKTEGER SGSA.SCjS.OAyr.FWGL 

CTeMF=(2»HAXT»-PTPIN)/3 

NTEMP* lPAxT+2«HI!<n/3 

IHLAl^LAl 

IF(LKLAI-GT.6.65J LHLAI=6,65 
IFtOAYT.EO.Q) GC TC S9 
LH=»1.0 

IFJDAYT.IT.SGSAJ GC TO 10 
IF (FLAGL.EO. U GC TO 12 
1F(LA I .LT.PLAM GO TO 13' 

GO TO 14 

13 LAIHA> = PU-I 
FLAGL=1 

12 LH=.75+.25*HAI/LA-IHAX3- 

14 IPIOA.rT.UT.SGSB J. GO TO 20 

Iim = . 1 065*SRy 177.08 + U69*{:HL All, 

GROSPF=2COO*( l-exP-(-.000CaS5* INTU.H'KS 
GO TO 30' 

20 INTL = .1065*SR»38.436.=»LFbAt**.4a6 

GROSPF=2000<C1-£XP(-.OCOC895*'1NTU )*-KS', 

GO, TO 30 

10 lHTL-=. 1065*SR*38.A36+L'FLAI^*.486' 

GRCSPH=2GOOFn-EXPl--OCOC59* INTL))*K‘S 
30 IFIHAXT.GT.33. I GC TC 16 
CO TO 17 

16 IFIHAXT.NE.HINT) GO TO 11 
TSTR=0-.0 

GO TO IS 

H TSTR = (PAXT-33)/.'( tPAXT-.S=*H INT- -S+PTH IN 1 1 

15 GR0SPH=GRCSPH*( l-TSTB) 

17 DR£SP=CL=t(.011667*CRCSPH+LH«CDH+t .002955+--.OOOl28*D.TEHPt.OOOQ67 

=* =fOTEPP*»2)l*(-i; 

NRE5P=(24-Ot.I*{,01 16 7*{G8QSPH+DRE’5P)+L:H*CCK* t.OQ2955+-.0001-23»- 
» NTEHP+.000067+NTEKF**2U*{-11- 

NCE=GRCSPH+-OR£SP+NReSP 
OH=,OC67^NC£ 

COK=CEJ'+CP 
99 PNTEHP=NT£PP ' 

RETURN 

ENO 

SUBROUTINE TCH { FCAGL ,L A I ,PLA I , K-, TACC ,6R0S PH,Sfi ,M AXT ,0R6SPi COM, 

♦ CL..NBESP.NCE.OH.HINT,PTNIN,LA'IHAX,XS J 

PEAL CAT .PCAi.TACCtGrCSPH.SR.MAXT.ORESP.CCMtXLf NR£SP,NCE.OH,NtNT» 

♦ CTEMP.NTEFP.CA IFAX , PTHI N , I LA [,XS 
REAL TC.TN.TSTR.TSTR1.TSTR2 
INTEGER FLAGl.K 

CTEHP = (2+HAXT + PTHtN + 3*273. 15) /3 
NTEHP=:(HAXT'+2»y'Ur + 3*2-73.l5J/3 
If IHAX1.GT.23. > GC. TC 20 
T0=23. 

TSIRl*^!. 

GO TO- 21 

20 TN=tPTFIN+MINT)/2 

IF (TN -GT, 23) IN = 23 
TD=HAXt 

IFITN.LT.5.) TN=5. 

TSTRl=l-(PAXT-23.) /IPAXT-TN) 

21 IFIPTPIN.LT.S.) GO TC 22 
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♦ F10.4.//.20X. <U (K)').... 

♦ FIO.4.//.20X.’ ALPH4 tP-T} •, 

♦ F10.4.//,20X» ‘ FALLRT - 

♦ FIO.4. //, 20X, 'SO tL CCHSTANT (HK DAY TO -t/2) 

♦ F10.4,//,20X,*X sue 5 (lUIT. MATER CG^Te^T IN 5 CM. LAYER, HT..'f 

4> FlC.4,/./.20X,'FIELC CAPACITY - 

♦ F10.4.//,20X, ' 6PTS FLLTIPLIER *, 

♦ Flo. 4) 

105 FORMAT! ‘O', 19X, 'PLANTING OATH... 

.20A4I 

106 FORMAT! 'C ' . I9X. ■ F 1 ELC/ LOCATI CN , 

*»..*,-20A4) 

107 FORMAT !I2,1X,I2,1X,I2) 

108 format ! 15F5.0) 

109 format (8E10.4J 

110 FORMAT!' ',33X,'TRAN POT. TOT. 1 NT. GROSS OAY NIGHT *, 

♦ •TOT. ‘,21X, 'CUH* ' RC CAY SR MXT HNT RAIN LAI EVAP EYAP', 
*• EVAP LIGHT PHCTC RESP R£SP RS5P NCE OM/TR OH •, 

♦ •OH 8PTS THETA KS'./,' (LYS) IC) ICl (PHI (HHP, 

♦ • tMHJ !PM) PE/CMCP (HG/HCH**2J7DAY}-I •, 

*• tHG/tCMt»2)) 

111 FORMAT!' I2,.I3,F6. C,F5.0,F5,0,F5. 1.F5.2 ,F6.2 ,F6 . 2 .F6 . 2,F7.0, 

♦ F7.1 .F7.1.F7-l,F7-l. F7.1.F6.1, F6.1 ,F6,l, F5.2,F7,l,F6.23 

112 FORMAT! 'OTOT.'.FI.O.VIX.FA.O.TX.FA.O.FS.O ,F5.0 ,F8.0 iFT.O, F7.0, 
♦F7.0, F7.0. F7.0.F6.C. F6.0 ) 

113 FORMAT!///, IX, laX, ‘TCTAL GROSS ', 36X, ' CUP' ,/, IX, lOX, 

♦•RAIN EVAP PHQTC TCT.RESP NCE OH/TRAN DM*, 

♦ /,1X, • HCNTH IMP) • , • ♦♦♦♦! PG/(T0H**2)/DAYn***^* , 

♦ 12X, • (PG/ICM^’fEJ ) • ,/ ,»+• .77! •_' 1 ,/-) 

114 FORMAT!' ' , 1 4, 5X, F5, 1. 4X ,F5. 1, 2X,F 7, 1 . 6X , F7. 1 ,.1X,F7.1 , 4X ,F5.l , . 

♦3X.F6.H 

115 FORHATI'O'.* T0T.',4X.F4.0,4X,F5.O,3X,F6,0,6X,F7.0,2X, 
♦F6.0.4X.F5.0) 

116 FORMAT!//,' • .IIX . • GftCSS > . 31X, • CUM TRAN TOT. RATIO*, 

♦ • RATIO RATIO FOT. ' , 15X, ■ INT. ' , / , • '.llX,' PHOTO TOT.*, 

♦ •RESP NCE DP/TRAN OH EVAP! 2) EVAP! 21 T3/EO!21 *, 

♦ •ET/EEI21 T3/T2I2) TRAN(2) SR LIGHT',/,' ',• 8HTS *, 

♦ •♦**!PG/! I0H»»21/CAY 1 ) ♦♦♦' ,9X, ' I H6/!CH*^21) ' ,51X,' I LYS) »> 

♦ •IPE/CPCP)* ,/, ' + < . 123! '_' ) ,/) 

117 FORMAT!' ' ,A4, A2 ,F 10 . 1 , F9-.1 , F9. I ,F9 . 1, F8, 1 ,F1 0.2-,F8.2 , F9.2 , 

♦ F9.2,F9.2,F8.2.F9.0, FIO.O) 

118 FORHATI-'O* , ' TCT. ' , FIO, 0 ,F9.0 , F9.0 , F9 .0 , 9X , F9 .0 ,F8.0, F9.0, 

*F9.0. F9.0,F8.0,FS.O.F1C.CJ 

119 FCRPATI//.' EMERGENCE CATE ! 8HTS=' , F4.2 , • OAYS= ' , 13, ' > . . ' » 1 2, 

♦ •/*,I2.'/',I2,/,* JOINTING OATS 1 BPTS=< , F4.2 , • OA YS=' , 13, , 
»I2.*/',I2,'/’,I2,/,' 3C0TING DATE ( 3HTS = • , F4, 2, * 0AYS = ',I3, 

♦ •1.,..',12,'/',I2,'/',I2,/,' HEACING OATS 1 8HTS= ' , F4.2 , • OAYS=', 

♦ 13 ,'l , 12. ■/' , 12,'/' ,12. /, ■ SOFT DOUGH DATE ! S MTS= ■ ,F4.2 , 

♦ •• 0AYS='.I3,' I.', !2,'/',I2.'/',I2,/,* RIPE DATE I BMTS= • ,F4, 2, 

♦ « 0AYS=',I3,') ',I2,»/' ,I2,'/',I2,/,'l'} 

120 FCBHATI813I2,F4.21) 

121 FORMAT !I3 ,1'X. 13 ) 

122 FORMAT!//, ' • , 19X , • I SRIGATI GN DATA ■ ,// / ,20X , ' DATE *,12X, 

♦ 'AHCUNT ICPl*,//) 

123 FORMAT!' ' ,19X, 12. '/' ,12.,"/' ,I2.16X,F5.2) 

1000 RETURN 

END 

SUBROUTINE S0RCPH1CAYT.HAXT,MINT.,PTP,IN.SGSA,SGS8, INTL.SR.LAl, 
•♦GRCSPH.ORESP.OL .CDH.NCE .NRESP ,0M , PLA I , LAI PAX ,F LAGL • PNT EMP, KS ) 

real paxi,hint.pipin,intl,sr,lai,grcsph,oresp,dl,coh,nce, 


OSS3IN4L Page p. 
^ POOR QUAim 



o o o on o o 


78 


TSTR2-K. 

GO TO 23 

22 TSTR2 = l-(5.,-PTf'l\)/(TO-PTHIN} 

TSTR=ISTR1*TSTR2 

23 ILAI=LiI 

IF IFLAGU-FO.l 1 GO TC 12 
IFIK-GE.3.ANO.LA I.LT.PLAT) GO TO 11 
GO TO 12 
U LAIKAX^PMI 
FLAGL=1 

12 GO TO (l«2,3,3.'«.4-].K 

4 IF (LAI. GT. 1,5) GO TC 13- 
RAI= aAl + .025 
IF (TACC.GT.4.S) GC TO 14 
GO TO 2 

13 ILAI = ILAI + ,025 

3 IF (LAI. L£. 1.5) GC TC 2 
1F( ILA I.GT.,4.62) ILA.I=4.62 

GRQSPH=( 1012 -101-2*6XP(-.00S5*SR*. 5739=' ILA 1**. 32961) *KS*TSTR 
GO. TO 17 
14’ GR0SPK=0.fl, 

GO T.O 16 

2 IF (HA XT.. LT, 5-0) GO, TO 16' 

GO TO 19 

18 GROSPH=0.a 
GO TO 17 

19 IFCILAI.GT.4.62) UA[=4.62 

GRCSPH = {'440-440#£XP(-.002S*SR*'.5739*ILAI*=»-32.96) J*KS 

♦ ♦TSTR 

17 IF (FLAGL.EO.l) GO TC 16 

OResP==(-.276Cl-.014775*GPQSPH-,002329i'CDH*(H-.0035087=«OTEHP+ 

* .0001ACTcXP*«2-))»OL 

NRESP= (24-OD* {-. 2 76C1- .014775=* { GRCSPH+CRESP )-.002329*CDH' 

♦ =*(l+.0035C£ *NTEPP + .0Q01’'N1GHP'»<=2) 1 
GO TO. 1 

16 DReSP=[-.27601-.Ol'4775*GROSPH-.002329*C0H*I ILA I/LA IHAX )* 

T (1 + .00350£7*CT EHP t.OOO l.=»0T£HP**2 )•) *CL 

NR£SP=(2.4-0L)=*l-.27601-.01'4-775*( GROSPHtORESP 1- .002329=*CDH* 

* ( ILAI/LA-IMAX)*(1+. CQ350a7=*NTEHP+.0001=)'NTEMP=!==i'2) ) 

.1 NCE=GRCSFh+CR£SP+hR£5P 

0H=NCE*.CX167 

CD('=CCK»CH 

RETURA 

END 

SUBROUT INE TRANS'! T I .HXH2C , TV, ICS . LA I , T , ALPHA ,RN , SSD i COVER.TAU* 
1 ALPHA V, DRY, SGYSCR, T2 ,TCCC, TACC, FALLRT) 

SUBROUTINE; CALCULATICN OF TRANSPIRATION- 


REAL LA 1 , HXH20 , XS . T2 ,TGOE ,-TACC , FALLRT 
INTEGER SCYSCR 
TAVAU = TI-(TV=»1S20.1 
KS=TAVAIL/ (FALLRT»HXF2C) 

IF (KS-GT-1..) NS=1 
IF IKS. LT. 0.0) KS=0.0 
GO TO (l.l,3,4),S0YSCR 

FCR SORGHUM ANC SOYBEAN 
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l TAU=E>P(-. J98*L‘AU 
IF'(LAIr2.5) 53.51.51 

FCR- VHEATr 

3' TAU=E)iP(--737*(.AIJ, 

IF (TACC.GT.3.0R.LAI-.GE.1.25) CG TO 30 
GO 10^ SO' 

30.- IF HACC.Gr.A'-) -GO TO 52 
GC TD 51. 


rCR' CORN. 


A DRY=t.'AI 

^TAU=EXP(-.3a9*0flYv0. 1A38J. 

IF, (I.AI.LE..38). TiU?l. 

IF n.AI.3c.2.3.GS.;7GCD.GE.lAOO.) CO. TO' Al- 
GC TO 50. 

Al. IF lLAI.Lr.2.5-.AAC.T6CC.GE'.2000) GO. TO 50 
GO TO-- 51, 


50' T=KS*At,PfAV‘'I1.0-TAUl*SSC»RN/58.3. 

GC’. TO 59. 

51 T=KS* (AL5H1-.TAU}*SSC»RA/58.3 
CO. TO 59 

52 T=KS*ALPr'AV»U-(-EXP( - .737*L.'AI )') J *S SC«RN/5 8.3- 
GO TO’- 59 

53' T=KS»AIFFAV*{ l-{ EXP (-.3S6^^LA-.in-),»SSD*RN/5 S.3, 

C. 

C 

99 IF fKS.eC'.O.;0)' GO TO S3 
T2=T*17KS 
GO, TO 9-7 - , 

53 T2=0.0 
97 .RE TURN 

END 

SUBPOUT I\c jmtL A'I.P LA I, INTL. GRUSPH , TRESP i.SCE .OH.CDH, SR ,ST ART, KSl 
REAL LAl.;SR,PLAt,IMl,GRCSPH,,TPSSp-,NCE,CMsCOH,KS' 

REAL' PAR'.PCTINT.ARG 
INTEGER STAST- 
IFILAT'.ET.O.O) GC' TO. 99 
START=-.l ' 

FAR=U-3»SS- 
IFILAI.GS.PLAI )• GO TCI 
ARG=10»LAI 

PC.TIMT=3.31*ALCG(ARG 1+63. 6' . 

GO TO' 10 

t IFaAl,.LE'.2.71 GC TG 2 
IFtLA I.Lr-5-2J GC TO- 3 
FCTINT=9S. 

GO TO 10- 

2 'PCT1NT=32..5*LA.I.+ .A6- - 
GO TG 10 

3 PCT1MT'=1 .2 7+L-AI+88 .A 
10 IN.TL.= PAR»FCTIMT*.or 

IFnUA'I.Cf.P'CAI.CR.LAI.GT.3.51 .OR.-INTL.GT'.SlSO.l GO TO 20 
GROSPH = ( 5-.5l£-.17+'[ATL»*5.056) *KS 
GO- TO 30 

20 GR0SPM = {,7l5+IliTL**,7Z8).tKS' 
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30 TRESP=.26*GR0SPH+A.455 

IFITRESP. LT.O.Oi TRe5P=0'.0 

NCe=GRCSPH-TReSP 

DH=.0C3*NCE 

CD«=CCP+0P' 

GO TO 100. 

99 IFISTABT-EO-OJ. GO TO ICO 
INTI_=C.O 
GR0SPl-=0.0 
TRESP=0-0- 
KCE=0.0. 

CHsO.O 
loo RETURN 
ENO 

SUBRCUTINE. POT EV A I LA 1 ,,CFLAG, RN, SR ALPHA , SSD , EQ , DA-YT J 
SUBROUTINE: CALCULATICN OF P,QT. EVAP., 


INTEGER CFLAG.OAYT- 
real LAI 
IF(CFLAG-3)5,6,7 
5 IF tLAI.LT,.3-0) GO, TO. 3 
IF ICFLAG-U ltl',2 

FOR SORGHUN 

1 RN = ..836a*SR-130, 78- 
go TO 4 

FOR' soybean 

2 RH=.aC49+SR-l35.-97 
GO TO 4 

"3 RN=.-7248*SR-5Q-U 
GO TO 4- 

FCR: wheat 

4 RN=,86-70*SR-163,56 

IF10AYT.LE.L68)RN=,9593*SR-213,10' 
iF(DAYT.GT.202> RH= . 925 8*SR-157. 4208 
GO TO 4- 

FCR CORN 

7 IF ILAI-GE.3-0J GC TC 8 
RN=.a6Ce*SR— 103.92. 

GO TC 4 

3 RN = .a4£*SR-U44.49- 

IF 1CAYT.GT.84) RN =. 766* SR-99.84 

4 EG'»ALPHA*SS0.*aN/58'.3 
RETURN 

ENC 


SUBROUTINE EVAPCFLAG.LAI ,RN , EST . £S , U.CNST .COUNT ,SS D< DRY ,C0VER» 
ICAYT, r AO.SOYSOR) 
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C SUBROUTINE; CALCULATICN OF SOIL EVAP. 

C 

C 

REAL LAI.RRNS 

INTEGER flAG.CAYT, SCTSCR' 

GO TO 12.6.6.5).SCYSCR 
C 

C FOR SCYSEAN' AND V.HEAr 

c 

6 IF (CCVER.GT.LAI) TAU=.eS2 
GO TO 2 
C 

C FOR CORN 

C 

5 TAO=EXP(-.38S*ORY+C. 14383 

IF (LAULE.O.SS’l-T-AU-t---' ' 

IF (OAYT. GT.90.ANC.LAI.lt. 3*67) TAU=.270 
C 

C SWITCH AS TC HHICH SOIL EVAPQRAT-ICN 'FORMULA 

C TO USE • 

C 

2 IF (FlAG-l) l,3',3 

I ES=TAL'*SS0+RN/58.3 
EST=£ST+es 

IF <EST.LE.U3 60' TC, 4 
FLAG=2 
£S =£5*0.6 
GO TO 4 

3 ES=CNST*(S0RT(CCUNT)'-SCRTtC0UNT-l.O) 1 
COUNT=CCUNT*l.O 

4 HRNS=T'AU*RN/S8.3 

IF IHF.SS.LT.SS I ES=HRNS 
: RETURN 

END 

C 

FUNCTION C£L.TA(T3 


CELTA=0.0155416*T - 0. OOOOOS*T**3 0.000 0001*7**4 0.40408273 

RETURN 

END, 


SU8RCUTINE CLCKER 

PURPOSE calculate PART OF BIO-TIME TODAYtBAIER MODEL 

OESCR-JPTKN -CF PARAH6IERS 

CCEF; COEFFICIENT TO CALCULATE TIMS 

TN: pin TEHP 

TY; hay TEHP 

DU CAY LENCTF- 

tacc; total of tihe farts 

TCAY: TCCAYS tihe PART 

Ns INCEX INTO CCEF 

SU8RCUT.INE CLOKERl COEF.-TN. TX .OL . TOAY .K . HULT . TACC ,C LOCK , HHO , J JJ t 

* YEAR, BTFLG.O IFF, SYR f SHC.SOAY.SCLK, ICOUNT, 

* RFLG,KAY,PTACC13 

c 

REAL CCEF <6, S3 , TN , TX ,CL , TDA Y ,PULT , T ACC , 0 I FF , KAY 16 1 .PTACCl 
INTEGER K, CLOCK. MHO, JJJ, VEAR,8TFLG,SYR( 63 ,SHO 16 !„ SO AY ( 5 3-,SCLK( 63 , 
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♦ ICCUNT.RFLG 

-c 

TX=M, 9 -A 5 .J*TX)+ 32 .- 
TN = ( C;./. 5 .)*TNJt 32 . 

C 

C FIND CAYtEhOTH CCNTB- 18 UT ION. TC .6 IC=-T XH 6 

C 

Vl=CnF(K,21’»(0L-C0EF(K,-n l.+CDEF IK . 3 )#{ OL-CQEFt K,.1J l**Z 
IFtVl.LT.C.) Vl=0-0 

Flf .0 HAX TEFP CCMRieuTICN TO. BIO-TilHE' 

V 2 =C 0 EF(K.Sli't TX-CCEF(K, 4).1 + CCEF( K, 6 )* (TX-C 0 EF.(K-, 4 ) i+=»a 
IF (V 2 .LT.O.O.CR,T.X-.IT. 23.641 V 2 = 0.0 

FIND FIN. TEMP CONTRIBUTION- TO BIQ-TIHE 

V 3 =CCtF{.K, 714 (TN-CC£FIK.„ 4 l ) <-COEF IK, 8 -)+ (TNj-COEF ( X , 4 n *+2 
IFIV 3 ..LT.. 0 . 0 ) V- 3 = 0 . 0 , 

TODAYS CCNTRIBUTTCN. TO BIQ-TIHE' 

TDAY=V 1 *{V 2 +V 3 I 

CUKHALATTVE sums for 6 HTS 

IF(K.£ 0 . 3 .GR.K.£e. 4 J TCAY=TDAY*«ULT- 
TACC=TACC+TCAY 
CLCCK=CLCGK+l 

IFUHFC‘.GT.. 3 -AND.HH 0 .LE,. 12 J .ANO.TACC.GT.I.BJ- TACC = 1,6 
■ IF.IBTFUG.EQ.U GC TC 10 

IFtK-£ 0 , 3 .AND.TACC.GE. 2 .TJ GC TO' 20 
GO TO 30 
20 . BTFLG = l 
GO TO 95 

30 CIFF=-TACe-PTACCl 

IF(S.YRt 2 )'.E 0 - 0 ./K 0 .TAGC.GE. 2 . 8 ), GO TC 40 
1 F( (TACC.GE.r.S 5 .AN 0 ,TACC-LT. 2.81 .AN 0 .X.NE* 3.1 GC TO 50 - 
GO TO 10 

50 IFt 0 IFF.G£-. 02 ). ICCUNT = IC 0 UNT.+ 1 
IF’( ICCUNT-UT.IOJ GO TO ICO, 

KAY.( 2 ),=TACC 
40 TACC= 2-00 
GC TO 9-5 

10 IFtRFLG.EC. IJ GC TC 111 

IFITACCIGE'-KAYIK-I.ANC.K. LT -61 GO TO 99 
IF IT-ACC.tT.. 5,001 GC TO 100 
RFLG=l 

99 SYR(K-)=YEAR 
SKClKjsHFC 
SOA,YIKJ=JJJ 
SCLKlKI=CLOCK 
CLCCK-C 

IFIK.LT .61 K=K'*l 
100 PT'ACC 1 ='TACC 
111 RETURN 
£N0 

c - 

c 

SUBROUTINE DAYUTHEVAL ,ZVAU,ORAIK,CO,TH£MAXl 
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SLBfiOllT.INe TC CCNTRCL, ORA INAGE- 


INTEGER CCJS) 

CIHENS IC.S THEVALCS) . ZVAL 1 5 1-. TAOO 15 J . THEM A)( ( 5 1 

CRA1N=0.0 

00 4 1=1,5 

TADC( n=c,o . 

TAOC t DiGS'A IN/ ZV ALM I 
TCK=ThEVAL( I-H-TACO(l) 

IF ITCK.L5..5) GC TO 1 
DRAINS (TClt-.S)*ZVAL l I) 

GO TC 2 
CHAIN=0.0 

IF (TFEVAU n.LE.TFERAX m J GO TO 4 

IF ICCU).LT.2> GO TC 3 

C0(II=1 

GO TO' 4 

CRAIN=(TFEVAL{ I I-THEFAXI I) )#ZVALM)+ORAIN 
-THEVAlU-JsTHEHAxm 
COUJsO 
CCNTINUE 
00 6 1*2.5 

IF ITACDI lI.EO.Ol GC TC '6 
CO ( I J=2 

THEVAL t UsThEVALU H-TACCfl) 

IF (TF£VALI'n.GT..5) ThE VALMI = . SO 

CCNTINUE 

RETURN 

ENO' 


SUBRCUTlXe CAYOI.TFSVAi:,RAIN,ZVAL,RUNCFF,CC,THEHAX> 
SU2R0UTINE TC CONTROL A RAIN., 


DIMESS-ICN ThEVAa5I,ZVAL.(5),.THEHAXI5) 
INTEGER CC(5) 

IFtRAIN.GT.25-.4l GO. TO 3 
RsRAIN- 
RUNCFF=O.C 
GC TO 4 

3 R=2S-4*{AAIN/2S.4I*».7S 
RUNOFF=RAIN-R 

4 CO, 5 1 = 1.2 

CK=t.5-TFsVAUI n*ZVALtn 
IF (R.LT.CK) go TO 6 
TH£VAUn = -5 
COM) =2 
R=R-CK 
5 CCNTINUE 

RUNCFF=R+?LNOFF 
GC TO 7 

6 THEVAL { IJsTFEVALI I ).+ {R/ZVALI D ) 

IF (TlEVALCU.GT.THcFAXt ID- COl 11 = 2 
7 RETURN 

ENO 

C 

SUBRCUTINE CISTRtT,TV4L,XVALl,KVAL2,LAI) 
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SU8R0UTIN£ TO CALCULATE THE TRAN SPI RAT ICN IN EACH- LAY.ER. 


CIHENSICK TVALI5) ,KVALUSi’,KVAL2l-S) 

SEAL KVAll .KVAL2 ,LAI 

If (LAI-GI.U GO TO 3 

CO. 2 1 = 1.5 

TVAL.1 1 I=T»KVALlt 13 

£0 TO 5 

CO 4 I?l,5 

TVALtn.=T»KVAt2(I) 

S RETURN 
ENO 

SUBROUTINE HOIST I THEVAL ,ES. TVAL .2 VAL .THE HINT 

TC CALCULATE THE SOIL MOISTURE CONTENT IN EACH LAYER 


01 HENS ICN THSVAL [ 5 ) -TVAL (5,3 , ZVAL ( 5 3 , THEM IN (5 3 

REAL ELEFTl.tLEFT2.RHC.ES 

£LEFTl-=eS 

IP (ELEFn.EO.0.03 GC TO 2, 

1 CO l<3 1=1,5 

IF- ITFEVALtn.EO..!) GC TO 19 
£LSFT2=cLErTl/,iVAU ( I 3 
7CK=TFevALni-£LEFT2 
IF (TCK.LT-.13- GO TO 11 
THEVAL (.tI = TCK 
GC TG 2 

11 £L£F-T2=TF£VALT I-3-. 1 

£LEFTl=eLEFTl-( ELEFT2*ZVAU 1 3 ) 

TK6VALtI3=-l 
19 CCNTTTiLE 


2, CO. 29 1 = 1.4 

3F tTVALM3.EO-0.0 1 GC TC 29 
IF (TFEVALMl.LE.ThEHINt II 1 GC. TC 21 
TCK=TFEVAL m-lTVALI 13/ZVALt 13) 

IF (TCK.LT.THEKIM n 3 GC TO 22 
THEVALt n=TCK 
GC TO 29 

21 TVALI 1 + U=TVaL( I + U * TVALU) 

TV AL( 11=0.0 

GO TC 29 

22 RMOsThEVAUn-THEHIHd) 

TVAL-1 IHT^TVALn + UMTVAU I )-i RHO=!‘ ZVAL (1 3 3 3 
TVAUI I )=RHC»ZVAL( M 
THEVAK n=-TFEHINn) 

29 CCNTINUe 


THEVALt51=TFEyALI5)-(TVAL( 5)/ZVALt-S) ) 

C 

RETURN 

ENO 

SUBRCUTINE CAY2t IFEVAL . I RRI G , ZVAL , RUNOFF , CO.THEMAX 3 
CIHENSTCN THEVAL(.5) . ZV AL( 5 1 . THEHAX 1 5 3, 
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REAL IRRIG. 
integer: CC(.5J. 

. r=iprig;- 

00, 1 1 = 1,, 3 

CK=( .S-Ii-e.VALtiU )*ZVAL'U J; 

IF^IR-. LT-VC'xi ■ GG JC 2 - 
THeVAEn‘J=^.5.- 
C0(:J) = 2 
R=R-CK 

1 ccntince: 

fiUNCfF=PtRUNOEE> 

GO. TC 3.- 

2- TH£VALM)=TE£VAD(-I )»R/2VAtU}' 

IF ( THEVAC t>U-,GT/.,TH£R AX-t I.) l-.CO fU=2; 
3 REJURN 
ENO, 
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